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Discrete Approximations to Elliptic Differential Equations’) 


By WoLFGANG Wasow, Los Angeles, California?) 


1. Introduction 


There exist many methods for the numerical solution by finite difference 
methods of DirRICHLET’s problem for linear second-order elliptic differential 
equations. Most of these use square grids and finite difference operators involv- 
ing four or eight neighboring points at each net point. But rectangular, trian- 
gular, and hexagonal grids have also been employed [1], [2]*). In more than 
two dimensions the number of grid systems that could be used, at least in 


' principle, is even larger. The emphasis on simple geometric patterns is natural 


from a computational standpoint, but this simplicity is by necessity destroyed 
near the boundary C of the given domain B, unless B is a very special type of 
polyhedron. The various techniques proposed for interpolation near the boun- 
dary all amount to a modification of the regular grid in some manner that 


depends on the given domain [1], [2], [3], (4). 


It is therefore desirable to develop the theory of such approximate solutions 


of DrricHLet’s problem in a way that is sufficiently general to contain most 
_ reasonable discrete schemes as special cases. The essential tool for such an 
~ investigation was developed by PETRowsky in {5}. In fact, PETRowsky’s 
~ approach is so general it includes also certain approximation methods which 


are not based on a fixed system of discrete points, such as statistical sampling 
(‘Monte Carlo’) procedures using random walks that are not bound to the points 
of any grid. 

PETROWSKY showed that his ‘generalized Dirichlet problem’ possesses 
a unique solution which approaches the solution of the ordinary Dirichlet 
problem, but he was not interested in the size of the error involved. Whe shall 


~ call this error—which is the generalization of what is commonly called ‘trun- 


cation error’—the discretization error. 
The first purpose of this paper is to study the discretization error in a manner 


- that sheds some light on the comparative value of the various approximation 


1) The preparation of this paper was sponsored (in part) by the Office of Naval Research, U.S.N. 
2) Institute for Numerical Analysis U.C.L.A., National Bureau of Standards. 
3) Numbers in brackets refer to References, page 96. 
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schemes which have been proposed. This will be done by a generalization of a 
method first used by GERSCHGORIN [6] and later modified by CoLLatz [7]. 

In section 5, these results are applied to an investigation of RICHARDSON’S 
scheme for the improvement of approximate solutions by means of a linear 
extrapolation from two results obtained with two different mesh lengths [8]. 
This method has proved its usefulness in numerous numerical applications [8], 
[9], [10], [11], but its validity in problems involving curved boundaries does not 
seem to have been investigated mathematically. It will be shown that, except 
for domains with very special boundaries, the benefit to be expected from the 
use of RICHARDSON’s idea in the numerical solution of DirIcHLET’s problem > 
depends decisively on the interpolation formula employed near the boundary C. 


2. The Integral Equation 


Let x, y, etc. denote points of n-dimensional Euclidean space E. Then all 
discrete approximations to DirICHLET’s problem considered in this paper are 
special cases of the integral equation problem 


sige - Uly)\ dP(yx), «eB, if 


f(x) ge gare. Cat 


Here and below the Stieltjes integral is formed with respect to the first variable 
in F(y, x), and B = B + C, C denoting the boundary of the bounded domain B. 
On Fy, x) the following fundamental condition is imposed: Be 

Condition I. F(x, y) is a cumulative distribution function in B, with respect 
to y. It is Lebesgue integrable in x and y combined for x € B, y € B. 

The function /(«) is supposed to be bounded and Borel measurable. 

This formulation is somewhat more special than that of PETROVSKY [5], 
who uses as kernels F(y, x) distribution functions over the whole space E. For 
practical purposes this means that our theory does not include discrete approxi- 
mation schemes to DirICHLET’s problem which involve points outside of B. 
Also excluded are random walk (Monte Carlo) schemes in which the moving | 
point can jump into E — B. A more general formulation including these possi- 
bilities would introduce no serious new difficulties, but as it would make the 
arguments somewhat involved, we shall limit the discussion to (2.1). 

The problem (2.1) contains as a special case all practically used finite 
difference approximations to elliptic differential equations of the form 


1 n 02 n 0 
Lul= 7 2, aul) ose + D, bit) = 0, (2.2) 
1,k=1 - j=1 I 


ern hs ane he 
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- where Xe (¢ = 1, 2, ..., m) are the coordinates of x. If a term in » itself is present, 


then the theory of this paper would have to be modified by dropping the 
condition 


, 


[ao At gS 


If the coefficient of wu is nowhere positive, this modification would be quite 
_ easy. It is, of course, possible to approximate equation (2. 2) by finite difference 


expressions whose coefficients vary in sign in such a way that condition I is not 
satisfied. But the convergence and stability of such approximations is doubtful, 
and they do not seem to have been used in practice. To any equation of the 
type (2. 2) a difference approximation satisfying condition I can be found. This 
result is contained in [12]. 

It should be noted that in the formulation (2. 1), the solution U(x) is defined 
at all points x, even if the Stieltjes integral reduces to a sum of terms correspond- 
ing to regularly spaced net points. For instance, if LAPLACE’s equation in two 
dimensions is approximated by the simplest finite difference expression the 
first equation in (2. 1) reads 


1 
U (ay, Xa) = 4 U(x, + hy, a) + Uy — he m9) + UH, a + Fh) 
(2.3) 
1 
Si ri U (x1, %, — hh), 


where x,, % are to be considered continuous variables. In other words, every 
point x = (%,, %2) generates here a square net of which it is a net point. 
In [5] PeTRovsky proved that problem (2. 1) possesses a unique bounded 
solution provided condition I and the hypothesis below are satisfied. 
Condition II. Let 


Fy, 2) =F, *), Ful ®) = [Fn aly) FG 2), m= 23,0, 
B 


then there exists a positive integer m and a positive number c <1, both independent 
of x, such that 


[aE nly: e Stn C8. (2.4) 
B 


If F(y, x) is interpreted as the transition probability of a random walk, then 
condition II means that a particle starting from x has a positive probability of 
leaving B on the m-th step (cf. [13)). 


i 
; 
m, ‘ 
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84 
In the study of the discretization error in the next section we shall need a 


few results concerning the operator 
£f] =/ u(y) dF(y, x) — v(x) . (2.5) 
B 


In writing x € B, etc., below, we always mean ‘for all x in B’. 
LEMMA 1. (Maximum Principle). If condition I is satisfied, then 


£[0] = | v(y) aF(y, x) — v(@) # 8>0, xEB,CB, 
B 
implies 
uP v(x) < sup v(x) . 
B-B, 


Proof. If the last inequality were not true we would have, for all x in Be 


Bee seP hel 1p a, x) + sup o(x) | dF(y, 2 sa — 6, 


‘teen eS 


which is*a contradiction. 
LEMMA 2. Let F(y, x) satisfy condition I, and let K, be a constant such that 


| Slol}=K,, forsen,c B (2.6) 
I} q(x) is a nonnegative function in B such that 
f[g]S—1, x€B,, (2.7) 
then 
sup | v(x)| < K, sup g(x) + sup | v(x)| . (2.8) 


B, B BB, 
Proof. Let K¥ > K,. From (2. 6) and (2. 7) it follows that 
Ela0 —K* 9) Ke KS 0 


and therefore, by Lemma 1, 


sup (a: — K¥ q) < sup (--v — K¥ q) S sup |v]. 


Pr BB, B-B, 


ie 


Hence, 
sup |v| Ssup |v| + K+ pe q, 


BeBe 


and (2.8) is obtained from this by letting K* tend to K,. 
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LEMMA 3. If in addition to the assumptions of Lemma 2 the three conditions 


[f[vJt[sk,, «€B—B,, (2.9) 

Wx = 0, x CC, (2.10) 

[aF,) Sk<1, we BoB, (2.11) 
B 


ave satisfied, then 


ju(x)| S —— (Ky, sup q(x) + K,). (2.12) 


Fa 


|o(x) | Ns / Iv \| dF(y,x)+K,, x€B—B,. 


- Depending on whether | v(x) | is larger in B, or in B — B,, application of (2. 10) 


Proof. The inequality (2. 9) implies that 
: 
; 
j and (2. 11) to this inequality leads to 


e sup |v(x)| < ksup |v(x)| + K, (2.13a) 
3 B=. B, 

. OF 

. sup |v(x)| S ksup |v(x*)| + K,. (2.13b) 
5 Bosh B-B, 


_ The inequality (2. 8) and, on the other hand, (2.13a) or (2. 13b) constitute in 
each case a pair of simultaneous inequalities for sup | v(x) | and sup | v(x) |, which 
B 


1 BB, 
- can easily be solved. We omit these short calculations, which in each case lead 
_ to inequalities that imply (2. 12). 


3. Bounds for the Discretization Error 


Let w be the solution of the problem 


hy es 02u 0 
Lul=z 3) 4) oy, d%, ey bg(r) 2% 0, eB, | 
1,k= z 
earners. et . | 


It is assumed that the domain B and the boundary function /(x) are such that 
this problem possesses a unique solution. The difference 


w=U—u, (3.2) 


86 WoLFcANG Wasow ZAMP 
where U is the solution of problem (2. 1), satisfies the integral equation problem 


w(2) = | w(y) dF (y, x) + g(x) , x€B; w(x)=0, x€C, (3.3) 


B 


g(x) = | u(y) aF(y, ») — u(x) (3.4) 
B 
From lemmas 2 and 3 of the preceding section we know that the discretization 
error w(x) is small, under suitable assumptions on F(y, x), provided g(x) is 
small. Our problem is then to find functions F(y, x) which make (3. 4) small. 
In all applications F(y, x) is considered to belong to a family F(y, x; pb) of 
kernels chosen in such a way that g(x) tends to zero as uw +> 0. 
It is quite easy to find functions F(y, x; uw) such that g(x) is as small as we 
please im any fixed closed subdomain B, of B. For instance, if L[u] = Au, we 
may take as F(y, ~; w) a kernel for which, in B,, 


[20) dF, 2»; w) 


s 


is the mean value of v(x) on the m-dimensional sphere of radius about x as 
center. Since u(x) is harmonic in B, the two terms in the right member of (3. 4) 
will then cancel for x € B,, provided yu has been taken small enough [cf. section 
4e)]. Hence g(x) = 0, for x €B,. This example illustrates the fact that the 
main difficulty of our problem resides in a narrow strip near the boundary 
C of B. 

Accordingly, we divide B into two disjoint subregions, B = B, + By, such 
that different properties of F(y, x; 4) are required, depending on whether 
x € By or x€ B,. The domain B, is bounded away from C. The assumptions 
which we introduce now are devised so as to apply at least to all those actually 
proposed finite difference techniques which use interpolation rather than extra- 
polation at the boundary. 

Let 


u(y) = P(y, x) + sr Bs (95 — 45) | u(x), x, yeB, (3.5) 


be TAyLor’s expansion, to terms of degreee  — 1, of u(y) about the point x. 
The function P(y, x) is a polynomial in ¥; (7 =1, 2,..., 2), of degree p — 1 at 
most, for all x € B. Also, 

P(R, Bie be). (3.6) 


The symbol % is to be treated as a variable in performing the indicated differen- 


tiations and then to be replaced by a certain point on the segment joining y 
and x. 


ee ey eR ee ee Se 


Vol. VI, 1955 Discrete Approximations to Elliptic Differential Equations 87 


Condition III. (a) There exists an integer p > 2 such that for the solution u(x) 
of problem (3.1) the corresponding function P(y, x) defined in (3.5) satisfies the 
integral equation 


P(x, x) - | Py, x) dF(y, x3) =0, x€B,. 
B 
(b) There exist constants M and «> 1 such that 


i) be p 
al (3 ls— il) dF (y, xp) SM yu", «xe B,. 
i \e 


Condition IV. For x € B, a statement similar to that in condition III is true 


| with the constants p, M, « replaced by q, N, B. Of q and B only the inequalities 
— g>0, B> 0 are required. 


In practical applications the kernel F'(y, x; ) is defined in such a way that 
the preceding conditions are satisfied for all solutions of L[w] = 0 defined in B. 


_ If, in particular, L[w] = Au the conditions are usually satisfied if P(y, x) is 


any harmonic polynomial of degree p — 1 in the y; (7 = 1, 2,..., ). 
The next condition is regrettably restrictive, but it is essential for our 


- method, and has been used in all previous work along these lines ({1], [2], [6], 


etc.). The condition can probably be relaxed, if entirely different methods are 
used, but the necessary arguments might be rather difficult. in [BEE e.g., [14] for 
a somewhat related problem.) 

Condition V. The p-th and the q-th derivatives of the solution u(x) of problem 
(3. 1) ave bounded in B. 

On the basis of these assumptions we immediately obtain bounds for the 
function g(x) of (3. 4). Let A, be a bound for the numerical values of all 7-th 


_ partial derivatives of u(x) in B. Then (3. 4), (3. 5), and (3. 6) imply 


Wane, ee Bus 
nea agen 3.7 


| 
sil ae 
( 


In order to apply lemma 3 to the integral equation (3. 3) for the discretiza- 
tion error, we need two more assumptions. 
Condition VI. 
[a re pyeh< 1, ve By, 
B 


where k is independent of mu. a 
Condition VII. There exists a function Q(x), nonnegative in B, such that 


£[Qj=-z, «eB. (3.8) 


Q(x) may depend on qu, but must be bounded as u > 0. 
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We shall not attempt here a general discussion of the existence and construc- 
tion of such functions Q(x). If the conditions of PETRowsky [5] are satisfied, 
which is the case in all applications known to the author, the existence of a 
nonnegative bounded function satisfying (3. 8) with the equality sign was 
proved in [15]. This function is the y-fold of the mean duration of a random 
walk in B with transition probability F(y, x; wu). If the operator 2 approximates 
LAPLACE’S operator, as 4 > 0, then, as was proved in [15], this particular choice 
of Q(x) satisfies the inequality 


[O)| Sat (PS 2) vA + elu)], 3.9) 


where V is the volume of B and the dimension. The function e(u) tends to 
zero with mw. An explicit construction of a function Q(x) satisfying (3. 8) and 
the other requirements of condition VII can be found in [1], provided 2(Q] 
belongs to a special but important class of difference operators in a square net 
in two dimensions. If this operator approximates LAPLAce’s operator, then the 
corresponding function Q(x) is numerically not greater than the square of the 
radius of the circle circumscribed about B. This bound is larger than the leading 
term of (3. 9), but it is valid for all yw. See also [19]. 

We now apply lemma 3 to the discretization error w(x). Thanks to (3. 3), 
(3. 4), and (3. 7) and the conditions imposed, all hypotheses of this lemma are 
satisfied. In particular, g(x) = w-! Q(x). Hence we have proved 

THEOREM 1. If conditions I to VII are satisfied, the discretization error w(x) 
of the generalized Dirichlet problem (2.1) satisfies the inequality 


)| Sy (M Aya sup Ole) +N Ay wt). 


ik 
\ B 


If no details are known about the relative sizes of M, N, A p Aq, it appears 
sensible, in most applications, to use functions F(y, x; w) for which a —1= 8. 
There is not much point in using, as is sometimes proposed, a very accurate 
approximation in B, unless it is matched by an equally accurate interpolation 
scheme near the boundary. 

The usefulness of theorem 1 is limited by the fact that u(x) is not known. 
This drawback is shared by the work of GERSCHGORIN and of CoLLatz of which 
it is a natural extension. In spite of interesting results by RosENBLOOM [16], 
Wats and Younc [17], [18], and oth< rs, no really practical appraisal in terms 
of the data alone seems to exist, except for very special regions. 


4. Examples 


It is not the purpose of this section to add to the approximation schemes 
suggested in the literature. The examples below are all quite wellknown. They 


Oa RENT Oe Se ne 


(U0) dF. 2) =F[U thm) +U—h) | 
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are inserted here in order to illustrate the preceding theory, and also as a prepa- 
ration for section 5. In all of them m = 2, and L(u) = Au. 
(a) GERSCHGORIN [6]. Let B, consist of all points x such that the four points 


(%1 + h, x2), (%1, %2 + h) are in B. For x € B, we define F(y, x) by 


: Rae TE 


+ U(4,, % + h) + U(x, x, — h)], | 


_ For x € B, = B — B, let &(x) be one of the points of C nearest to x and define 


[ U0) aE, ) =U), xe Be: 


= 


_ then assumptions I to VII are satisfied, with 


ek ht op 2AM 


al 
ge = Ag tye NS Pay ke 0. 


The inequality of theorem 1 becomes 


\w| red oe x) +Ayh. 


(6) CoLtLatz [7]. The only change from GERSCHGORIN’S scheme concerns the 
definition of F(y, x) in B,. If « € B, then at least one of its four ‘neighbors’ 
(x, + hy %2), (%1, %_ + h) lies in E — B. On the segment joining one such point 
to x let x1 = x1(x) be a point of C having a distance h,<h from x. If this seg- 
ment is extended by a distance h, < h beyond x a point x* is reached, which is 
either the first intersection with C on this side or the nearest net point of the 
lattice, whichever comes first. Then F(y, x) is defined in B, by 

JG) aly, a= Be By (4.2) 
B 


In other words, we use linear interpolation between x, x! and x*. Here we have 
1 
g2 ed i=l, p= 1, k=. 
Hence 


jw] < (gy damax QC) + 2 As) he 
B 


(c) Mikeladze [3]. Again only the definition of F(y, x) in B, is modified. 
Consider, as before, the four lattice ‘neighbors’ of x. A neighbor which lies 
outside is to be replaced by a point of C on the segment joining it to «. Thus 
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four points x (7 = 1, 2, 3, 4) with respective distances h; from x are obtained. 
Let the numbering be done counterclockwise. Then we define 


[%) aL (Mk) = a ais igh me ey 
B 


U(#*) U(#*) U (#4) 
Nig(hg+ hy) — Iig(g+ hg) — Ag(hg + My) I? 


x€B,, (4.3) 


and find 


3 3 1 2 
q= 3, Ne =o aes se a ea 


(d) BicKLEy [20], MIKELADZE [21], MILNE [2], etc. The ‘weights’ defining 
the step function f(y, x) in B, in the previous examples can be illustrated, in 
a self-explanatory manner, by the symbolic pattern 


ae 
Rites 


@Q“@ 


1 


4 


(sometimes called a ‘stencil’). The stencil 


defines, in analogous fashion, a function F(y, x; 4) in B, for which 
p=, p=8, 2-4, M=0-000622. 


The definition of B, has to be changed now so that it contains only those 


points of B for which all eight neighbors of the stencil above lie in B. The use 
of this complicated F(y, x) in B, is justified only, if it is coupled with a defi- 
nition of F(y, x) in B, which contributes a similarly small error term. Such 


definitions can be given in many ways, but they always involve considerable 
numerical work (see [2], page 150). 


ea eee Ny a ee ee eee ee 
J tw ” . » Eke 
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where K,, is a circle of radius about x. This choice of the function F(y, x; 1) 


was ment:oned in section 2. Here B, is to be defined as the subdomain of B 
whose points have distance greater then wu from C. The first right-hand term 
in (2. 21) can here be replaced by zero. But since any definition of F(y, x; ) 
in B, which gives good accuracy near the boundary is necessarily complicated, 
the advantages of this kernel for Monte Carlo or other numerical procedures 
are probably illusory. 

5. L. F. Richardson’s Extrapolation Procedure 


Let U(x, w) be an approximate solution of a problem whose true solution 
is u(x) = U(x, 0). If U(x, p) is differentiable with respect to uw at w= 0, Le., if 


U(x, uw) = u(x) + (x) w+ off) , (5.1) 
then a knowledge of two approximations 
U(x, w,), U(%, wa) with w= 04, O<0<1 (5.2) 


can be used to construct the approximation 


Ee eI ls) + of) , (5.3) 


which is an improvement, as far as the order of magnitude of the error is 
concerned. 

RICHARDSON, in [8], recommended this procedure for use in connection with 
the solution of differential equations by finite difference methods. It does, 
indeed, lead frequently to striking improvements of the numerical results (see 
[8] and [11]). All these successful applications of the extrapolation are, however, 
concerned with problems in domains so simple that the problem of interpolation 
of the boundary values does not occur. The domains in question are almost 
always intervals or rectangles. 

In order to understand the decisive role of the boundary in this problem let 
us analyze in detail two extremely simple illustrative examples. 

Example 1. Consider the technique described in section 4, example (a), in 
the one-dimensional case. Let B be the interval | «|< 1, and let the prescribed 
boundary values be f(1) = 1, f(—1) = —1. The differential equation problem 


uw" = 0, |x|<1, u(—1)=—1, wfl)=1 (5.4) 
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has the obvious solution u(x) = +. The corresponding finite difference problem 
is here 


2 


Os) 1 Lh eS Ue ee he ee Te 


1 U(x +h) +> Ulx—h)—U(x)=0, |x| S1-h, | Be 


This formulation is the one used in our preceding general theory, but it is 
more complicated than is necessary in this simple example. The natural proce- 
dure here is to consider, from the outset, only values of h, for which 1 /A is an 
integer. In this case x could be limited to points of grids symmetric to the 
origin and having net points at + 1. But let us assume, for purposes of illustra- 
tion, that the grid chosen, although symmetric to x = 0, has a net point at 
1 — 6, with 0<6<h. Then the corresponding solution of (5. 5) has the value 


U(x) = U(x, h) = 57-5 =a t 2 wht O(n) (5.6) 
at the points of this grid. 

Now apply the extrapolation technique to this solution, setting successively 
h = h, and h= hg = 0 d, with corresponding values 6 = 6,,°0 = df, li h, and 
hy are chosen so that the same point x is a net point of both grids, formula 
(5. 6) is valid in either case, and we obtain 


U(#, hy) — 0 U(x, hy) © Soo Gat 


= ae ee hy thy PNB) 

Thus, the order of magnitude of the error is not changed by the extrapolation. 

In fact, for many choices of h, and h, the error will even have been increased. 
We mention, incidentally, that the explicit solution of problem (5. 5), fOr 

all x and his 


nit at 5 OL aie Tt = xh 
Lila ET ty fe 1+ #/h] + [1 — #/h]’ 


where [z] denotes the integral part of z. This expression can be used for a more 
detailed study of the error term. 

In domains in more than one dimension with curved or irregular polygonal 
boundaries some interpolation of the boundary values is unavoidable, and the 
preceding example forces us to the conclusion that there is no benefit. in the 
extrapolation technique, if the treatment of the boundary is no better than in example 
(a) of section 4. 

The same statement can be made with regard to the improved boundary treatment 
im example (b), section 4. This cannot be illustrated by the example (5. 4), 
since the discretization error is strictly zero in this case. Since two-dimensional 
counter-examples would be too involved for explicit treatment, we choose again 


. Vol. VI, 1955 Discrete Approximations to Elliptic Differential Equations 93 
4 


2 


ia leh 


la iS A ee eek 


a one-dimensional problem, which, although it is no longer a special case of 
_ problem: (3. 1), illustrates well the situation that must be expected in general. 
The boundary difficulties are the same as those encountered with problem 
~ (3. 1) in more than one dimension. 

Example 2. We set, for abbreviation 


sinhl=a, coshl=d. 
Together with the continuous problem 


“"—u=0,-|x| <1, Whee, u(—1) =— —, So) 


we consider the discrete problem 


BU Ph) + > Ue) UG) — 4 U0, |x| s1—-h, (6.88) 

a SR a ae eT es ee ae 5.8b) 

ey ies oe Sy ae soe! 2 i 
hy ' h 1 A 

pig Ulta yg U)=0, -1Sa<-1th. (6.80 


Bxere 1 — h, and A, — 1 are the right and left end points of the largest interval 


that can be fitted into |x| <1 so that x as well as those end points are grid 
4 points of a net of mesh length h. 


In analogy with example 1 we restrict ourselves to a special grid, symmetric 


f to the origin, with a net point at x = 1 — 6, 0<h. Because of the antisym- 
_ metry of the discrete problem and the symmetry of this grid, 


Ula) = -U(-2) 


- for points x of the grid. All solution of (5. 8a) with this property are of the form 


U(x) =csinhpx, |x|S1-—4, (5.9) 
where 
CAT ee Gas (5.10) 


2 


We now determine the constant c in (5. 9) in such a way that (5. 8b) is satisfied 
for the grid points in 1 —A< x <1, ie., for the one point 1 — 0. For reasons 
_ of symmetry (5. 8c) is then automatically satisfied. With the notation 


re) 
y= h 


the value of c turns out to be 


es tt 4 y) sinhp (L— yh) —ysinhp[1—(y+1A]}*. (5.11) 
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As (5. 10) implies 
h2 
p =1— eng O(h4) , (5.12) 


we find from (5. 11) and (5. 12), after a short calculation, 


1 
c=1+ [Ze + Fy (t+y)| A+ 04). (5.13) 
Insertion of (5. 12) and (5. 13) into (5. 9) followed by an expansion in powers 
of h? leads to 


U(x) = sinh x + [A(«) + B(x, h)] h? + O(h4) , (5.14) 
where 


: b 1 : 
A(x) = 7 (sinh x —coshx), B(x, h) = (2. ca a ed 7) sinh x . 


The first right-hand term in (5. 14) is, of course, the solution of problem (5. 7). 
If x + 0, the quantity B(x, h) in the first error term is a discontinuous function 
of h without limit as h + 0 because of the presence of the ‘sawtooth’ expression 
y. It follows, exactly as in example 1, that no improvement can be expected 
from RICHARDSON’S extrapolation. 

If the interpolation at the boundary is made more accurate than CoLLATz’s 
method used in the preceding example, an affirmative statement can be made 
concerning RICHARDSON’S extrapolation scheme. 

THEOREM. Let B be a two-dimensional bounded domain with boundary C, and 
let u(x), u4(x) be the solutions of the problem 


Au=O0inB, u=fonC, 
(5.15) 


Ome Ot Ne. 
Au, = ees ae 4 in: B,~«“=0OonC, 
respectively. Assume that u and u, possess bounded sixth and fourth derivatives, 


respectively, in B. Denote by U(x) the solution of the finite difference problem 
(2. 1), with F(y, x) defined by (4.1) in B, and by (4. 3) in By. Then 


U(x) = u(x) + uy(x) h® + w,(x, h) h3, (5. 16) 


where w,(x, h) is uniformly bounded in B, as h +0. 
Proof. Insertion of (5..16) into (2.1) leads to 


w(x, h) = [e(y, h) dF (y, x) + G(x, h) , (5.17) 


B 
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ra 
ae 


)=h-* | July) dB, 2) = uta] 


B 


(5.18) 


se | | m(y) aF(y, ») - a] =h-® Q[u] + ho? Of]. 
B 


With our choice of F(y, x) this becomes, for x € B,, 


[VAC 04u h 
G(x, h) == tee A eae! 7 Au, + 0(4) = 009), xEB,, (5.19) 


because of (5. 15). 
For x € B, it follows from (4. 3) and Au = 0 that 


L[ujl= 0), «eB. 


Ry RE ee NE 
i?) 
rae 
s 
> 


eer ee 


We also see directly from (4. 3) that 

: |u| = O(h?), «EB. 
: Insertion into (5. 18) produces 
| Giz, h) = Ol). ¥E By. 5,20) 


_ Thanks to (5. 19) and (5. 20) and the properties of F(y, x) stated in section 4(c) 
the function w,(x, 2) satisfies the assumptions of Lemma 3 in section 2, with 
_ K, and K, independent of h. It follows that w,(%, h) is indeed bounded as h + 0. 
| The preceding proof shows also that the largest contribution to the third 

right-hand term of (5. 16) — at least in order of magnitude with respect to h— 
originates in the boundary strip By. 

It is clear that the technique used in the proof of theorem 2 could be applied 
- to numerous problems: other differential equations and other discrete approxi- 
- mation schemes could be treated; if u(x) has sufficiently many derivatives in 

B and if F(y, x; «) depends sufficiently smoothly on mu in By, the formula 
_ (5. 16) could be replaced by an expansion to more terms. The author doubts, 
however, the numerical usefulness of such investigations. 

Summarizing the meaning of theorem 2, we can state that the order of the 
error term in the solution of Drr1cHLeT’s problem can be increased from 0(/*) 
to 0(h3) by an extrapolation using two values of h, provided the boundary 
interpolation described in 4(c) is used. Since this interpolation formula is fairly 
complicated, it is quite possible that one calculation using CoLLATz’s scheme 4(5) 
and a smaller value of / will yield the same accuracy with less labor. This 
criticism does not apply to applications of RICHARDSON’S extrapolation to 


problems on line segment or rectangles. 
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An extension of the methods of this paper to boundary-value problems of — 
second and third kind would be of interest. Some of the results of BATSCHE- _ 


LET in [22] go in this direction. 
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Zusammenfassung 


+ 
J 
; Der Fehler, welcher entsteht, wenn man das erste Randwertproblem fur 
_lineare elliptische Differentialgleichungen zweiter Ordnung mit endlichen Diffe- 
-renzenmethoden lost, wird mittels einer Methode abgeschatzt, die so allgemein 
ist, dass sie fast alle in der Praxis auftretenden Netze und Interpolationstechniken 
an gekriimmten Randern umfasst. Diese Abschatzung wird dann bentitzt, um 


die von L. F. RicHARDSON empfohlene Methode zur Verbesserung der Nahe- 
-rungslésung durch Extrapolation mathematisch zu analysieren. 


: (Received: April 22, 1954.) 
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On the Stability 


eet ae Oe 


j of the Equilibrium of a Linear Mechanical System 
By OENE Bottema, Delft’) 


: 


| In some papers H. ZIEGLER?) discussed the stability criteria for mechanical 
systems. In particular he drew attention to the fact that for non-conservative 
_ systems one has to be very careful with the application of other criteria than 
the kinetic one. Moreover dealing with such systems he encountered some phe- 
nomena which at first sight show a more or less paradox character. We mention 
for instance his discussion of a certain non-conservative system with two 
degrees of freedom, for which the condition of stability changes in a disconti- 
-nuous way when we pass from the case in which the friction is very small to 
‘that where it is wholly neglected. In what follows we add some remarks 
to the problems in question. We consider the general linear dynamical system 
_ with two degrees of freedom and discuss stability under the conditions just men- 
tioned, meeting 7. a. ZIEGLER’S result and answering the question when the 
said discontinuity exists. Moreover we state another theorem in the same field, 
which shows that the collaboration of a non-conservative force (a circulatory 


1) Technische Hogeschool. 
| 2) H. ZrecLER, Stabilitdtsproble 
Die Stabilitdtskriterien der Elastomechamik, Ingenieur 


bility, ZAMP, 4, 1-50 (1953). 


me bei geraden Staben und Wellen, ZAMP, 2, 265-289 (1951); 
archiv 20, 49-56 (1952); Linear Elastic Sta- 
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force in ZIEGLER’S terminology) and a gyroscopic force is fatal for stability | 
an equilibrium however strong it seems is ruined by such a combination 
however small the participants may be. At last we prove that the same theorem 
holds for a system with three degrees of freedom but for certain exceptions. 
.) 


2. 


We consider a holonomic scleronomic linear system with two degrees of 
freedom, the coordinates x and y of which are chosen in such a way (which is 
always possible) that the kinetic energy is 


Tas (e+ 9%). ) (1) 


We confine ourselves to small displacements from the equilibrium configuration: 
x = y = 0. Hence the Lagrange equations are as follows 


t 


K+ 4y,% + a,yY+dy,%+ dy =0, | 


: . (2) 
Y + Ay X + Ag ¥ + bo, % + bog y = 0, | . 
where a;; and },; are constants. | a; || is the matrix of the forces depending ont 
the coordinates, || };;|| that of those depending on the velocities. If | a;;|| is: 
symmetrical, there exists a potential energy function “y 


il 
Visese (ay, %? + 2 ay, XY + Ags y?) 


if it is antisymmetrical, the forces are circulatory. When the matrix || b;3|| is: 


symmetrical, we have a non-gyroscopic damping force, which is positive when. 
the dissipative function . 


if . caer “5 
2 (by, 4° + 2 dig % y + dos y*) 


is positive definite. If | },,| is antisymmetrical the forces depending on the: 
velocities are purely gyroscopic. 


| 4@,;| and || 6,,|| can both be written uniquely as the sum of symmetrical and| 
antisymmetrical parts: 


where 
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Bvhere 
it al 
bs =b4,= a (O15 + by), O= 2 (D1. — by). 


FR RRR ee 


We can simplify our equations by a suitable transformation 
%=%,cosy— y,sing, y= x,SINY+t 9, Cosy, 


which does not alter (1), leaves the type of equations (2) the same (only trans- 
forming the coefficients) and for which ~ and @ are invariant; we can choose 
g in such a way that either a¥, or bt, vanishes, reducing the potential energy 
or the dissipation function to a sum of squares. 

We take the first of these possibilities and writing again x, y, by,,... for 
1, V1, b%, etc. we get the following set of equations of motion: 


#+ayxtpyt bux t (be +o) y=, | 3) 
ee ey yet. 2) tes = 0_| 


The system has a potential energy function when # = 0, it is purely circulatory 
for a4, = 4%, = 0, it is non-gyroscopic for w= 0, and has no damping when 
bi = O42 = bo = 0; if damping exists, we always suppose that it is positive. 


eT ee ee oe To See arn o 


3. 


s'il ellis 


In order to solve the equations (3) we put x= C,e"", y= Cy e** and we 
obtain the following equation for the frequencies 7s 


A274 by,A+ ay, (by + 0)A+ 


| | =0 (4) 
4 | (dy, — @) A— AP + Boy A + Ag. 
or ' 
Coat 4 cy A® + 6, A+ CA + Oy = 9, 
where 

C5 = 1, 

Cy = by + Doo, 

Cy = Ay, + Ag9 1 (O11 Or, — Bia) a ee > (5) 


Cy = Ayy bag + Aon dy, + 2PO, 
2 
Cy = Ay, Ag, + P°- 


For ‘the equilibrium to be stable all roots of (4) must have real parts which 
are <0. For the latter condition it is necessary that c; 20 (¢ = 1, 2, 3, 4). 
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If c; > 0, we have moreover the necessary and sufficient Routh-Hurwitz 
condition , _ 

Co CR + C4 C3 
If cy = c,; = 0 the equation has two pairs of opposite roots and the condition 


reads 
CF 4 Cy Can (7) 


If, however, c, = 0, cz + 0, the sum of the roots is zero ; it is therefore obvious 
that they cannot have all a non-positive real part unless they are all purely 
imaginary, but that is impossible in view of Cz + 0. 

That means: for c, = 0, cs; + 0 the equilibrium is always unstable. (The 
same holds of course for c, + 0, cz = 0). 


We consider some special cases. 
(i) Suppose that the field force is conservative (6 = 0) and that there is no 
damping (b,, = by. = by, = 0). Then 


cy = 0, Co = @y, + gg + w?, c,=0, Cy = Ayy Ago 


so that we must make use of (7). If for @ = 0 the equilibrium is stable (a,, > 0, 
4, > 0) it remains stable when gyroscopic forces are added. If it has one 
unstable degree of freedom (a,, < 0, Ay, > 0) it remains unstable under gyro- 
scopic forces, because c, < 0. If both degrees of freedom are unstable (a; <0; 
@y_ < 0), the system can always be stabilized by a gyroscopic force, the con- 
dition being | w| > V— a,, + V— fee 

All these results are special cases of wellknown theorems, concerning the 
stability of conservative systems with an arbitrary number of degrees of 
freedom!), 

(ii) If the field force is conservative (p = 0) and if there are no gyroscopic 
forces (w = 0), then we have in the case of damping 


oa ries 2 Re. 
Cy = by + doo, Cy = Ay, + Age + (Dy, dag — 02,), C3 = Ay, Dog + Ago byy, C4 = Ay] Ago. 


Hence c, + 0, cs + 0 and SUPPOSING 44, bo. + agg by, > 0 the inequality (6) 
reads: 


(441 Dep + ayy dy1) (P11 + Dos) (B44 bee — dix) + by bys (441 — gg)? > 0. 

1) W. Tuomson and P. G. Tait, A Treatise on Natural Philosophy, Vol. 1, Part I (Cambridge, 
1879), p. 370-415, — H. J. E. Beru, On the Stabilization of Instable Equilibrium by Means of Gyro- 
scopic Forces, Phil. Mag. 49, 447-450 (1925). 
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That means that if the equilibrium is originally stable (a,, > 0, ag. > 0), it 
‘emains stable when damping is present. If it is unstable, it cannot be stabilized 
by the damping forces; for if a, <0, dg, > 0, we have c, < 0, and if a,, < 0, 
i>. < 0, we have c,< 0. 

_ (iii) If we have the same condition # = 0 as in (ii), but gyroscopic forces as 
well (w + 0) the coefficients (5) do not change but for c, which increases. 
Hence stability or unstability are not affected by the gyroscopic force. If the 
‘system is originally stable (a,, > 0, a. > 0) the situation is not changed by 
gyroscopic and damping forces; if it is originally unstable, it may be stabilized 
by a gyroscopic force as seen in (i), but if a damping force is added, however 
small it may be, the equilibrium is again unstable’). 

(iv) We now suppose that the field forces are non-conservative (p + 0). If 
damping is absent (b;; = 0), then 


= = 2 e = 2 
4 Cy = 0, Co= Ay t+ Ag+ Ww, Co= 2PM, Cy = Ay Ag2 + P*. 


. We have seen above that for c, = 0, cs + 0 the equilibrium can never be stable. 
Therefore: In a frictionless system there can be no stability if a non-conservative 
"force and a gyroscopic force (however small they may be) are acting together®). 
_ _(v) Suppose that we have a non-conservative system (p + 0), with damping, 
_ but without gyroscopic forces. Then: 


Suppose now that the damping force decreases in a uniform way, that means 
take b,,=€0,,, big = Dig, Doo = € bin, where bj,, D2, by, are constants and 
_e-> 0. Then for the inequality (8) we get 


4 Cy = Dyy + Don, Cy = Ayy + Aaa + (011 Dog — Oia), 

“ 

4 

q Cy = Ayy Dyy + Mag Dy, Ca = G11 Mog + DP”. 

Assuming that a,, > 0, ay. > 0 the condition for stability reads 

2 b bis + bag)? 

4 2, (un Doo + eae by3) (Ox boo — bis) (b11 + 529) 2 (Ou + Pas) i 8 
; (44, — 9)” + by, Oo0 wai by1 Ogs (8) 
| 

j 


‘ bf bas) 
(441 — pp)” > g, with 4 = se ae D7. (9) 


11 Obs 
- But if there is no damping, we have to make use of condition (7), which gives 


. (a41 — 40)” > &o where gy = 4 p?. (10) 


1) W. THomson and P. G. Tarr, l.c., p. 391: ‘Gyroscopic forces cannot convert instability 


; into stability if there is. any dissipativity’. 
| 2) If m and p are both small, the real parts of the roots are + w p/ (ay, — 499) for a4, + Ag9, 


and + p/(2 Vay) for a1 = 49; in all cases there are roots with a positive real part. 
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Obviously 
: (511 — 592)? = 
£1 — £0 = = (). 


by 22 


That means that in all cases where b,, + by, we have a discontinuity in our 
stability condition, such that the equilibrium may be stable if there is no 
damping, but unstable if there is damping, however small it may be. 

This is the result which ZIEGLER found in a special case. 

(vi) If the field is non-conservative and there is a gyroscopic force and a 
damping force as well, we have the general case. It is to be expected that here 
stable and unstable systems exist depending on the values of the parameters. 
This is confirmed by an example. Take 


Gy, =44=—4 > 0, = bes = 6, by. = 0. 
Then we have 


G=2b, &=2a+b*+m*, ce=2(ab+po), q=a?+ p?. 


Stability occurs when 


2 2 
y= adit @ > 0. and? ep eet eee 


C1 Cg 
which taken together gives rise to the single condition 


ab*+pmb—p2>0. 


As we have seen above there is no stability for b = 0; for small values of } 
there is no stability either, but if } is sufficiently large, the equilibrium is stable. 
The condition is 


2ab>—pw+Vp%(o2+ 4a), 


which is remarkable in so far that it does not remain the same if we change 
one of the parameters # and w into — p or — w respectively. 


5. 


_ In (iv) we stated that for a frictionless system of two degrees of freedom no 
stability is possible when there is a circulatory force and a gyroscopic force as 
well. We now ask whether this theorem also holds for a system with three 
degrees of freedom. If we take the kinetic energy and the potential energy of 
the conservative force as a sum of squares, then, using analogous notations as 


a a 


y 


a 


NN yee Ney 


«+ 


oe ses ee 


a 


” 


xe 
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above, we get the following frequency equation 


| AA + dy A + a, (Big + O3)A+ fs (big — Wa) A — Po 
| (bo, — @3)A— py A® + DogA + ay (bop + Ws) A+ Py| = 0. (11) 
| (bg, + We) A+ po (032 — 1) A — py A? + O534 4+ ag 
Here a, are the components of the conservative force, p;, w; and b,; denote the 
circulatory, the gyroscopic and the damping forces respectively. If the latter 
is absent the equation reads . 


dy A® + d, dA®° + dy At + dy A434 dl? + ds44+ d,=0, 
where 
a,=1, 
d,= 9, 
dy = 1 + dy + Gy + wi + w3 + @5, 
dy = 2 (py ©, + pz 2 + Ps Ms); 
dy = PRA PRA PR + ay dg + Ag Gy + Ay Ag + Ay W} + Gy WR + Ay 3, 
ds = 2 (ay py ©, + ae py We + 43 D3 Ws), 
dg = Ay Ay Az + A, Pi + Ay PE + Az 5. 


- The sum of the roots is zero; therefore, if they all have non-positive real parts, 


_ hence 


only purely imaginary roots exist. This implies 
dA? +d;=0 and A®°+d,A4+4,/?+4,=9, 


which is not an identity. Therefore the systems are in general not stable. There 
may be stable ones, but then the parameters must fulfil an equation (and not 


| only inequalities). It is not difficult to construct an example where stability 


is ascertained. Take 


ee en 4 ay 22 ag 2, Pa 4) OL Oy Os 2, 


then 
oo 24, d,=0, d,=—116, ¢,=0, d= 144 


and the equation is 
(32 + 2) (22 + 4) (A? + 18) = 0, 


which has no roots with a positive real part. 
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Zusammenfassung 


Anschliessend an die Zieglerschen Arbeiten wird das allgemeine lineare dyna- 
mische System mit zwei Freiheitsgraden in bezug auf seine Stabilitat untersucht. 
Es werden dabei nichtkonservative Feldkrafte, lineare Dampfung und gyro- 
skopische Krafte in Betracht gezogen. Die Anwendung des dynamischen Stabili- 
tatskriteriums gibt einige Satze, welche iiber die schon von THomson und Tait 
herriihrenden Resultate hinausgehen. Es zeigt sich, dass in einem nichtkonserva- 
tiven System ohne Dampfung bei Anwesenheit gyroskopischer Krafte das Gleich- 
gewicht niemals stabil sein kann. Dieser Satz gilt fiir Systeme mit drei Freiheits- 
graden nicht ohne Ausnahmen. Ein von ZIEGLER gefundenes Ergebnis, wonach 
eine Diskontinuitaét in der Stabilitatsbedingung auftreten kann, wird allgemein 
diskutiert. : 


(Received: February 24, 1954.) 


An Electromagnetic Difference-Equation of Importance 
in the Theory of Traveling-Wave Tubes 


By HErnricH DERFLER, Ziirich!) 


§ 1. Introduction 


It was shown [1, p. 23]?) that the currents and voltages on a chain of iden- 
tical three-terminal-pair junctions (Figure 1) are interconnected by the follo- 
wing difference equation: 


Vin +1) — 2ch IV (n) + V(n —1) = Zsh I'I(n) 


—W{[I(n +1) —2ch PI (n) + I(n —1)]. . 
The parameters J’, Z, W can be expressed in terms of the impedance para- 
meters of these junctions. 

It will be shown in this paper that 
equation (1) holds even under more 
general conditions, including inter- 
actions over long distances such as 
they occur in traveling-wave-tube 
structures, provided the parameters 


: Figure 1 
I’, Z, W are suitably redefined. Chain of three-terminal-pair junctions. 


1) Institut fiir Hochfrequenztechnik der ETH. Paper supported by a research fellowship of 
the Battelle Memorial Institute. 
2) Numbers in brackets refer to References, page 114. 
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§ 2. The Representation of the Traveling-Wave Interaction 
by an Integral Equation 


The motion of a cylindrical electron beam under small signal conditions 
- leads, by the Maxwell-Lorentz-field theory, to the following integral-equation. 


1 7 
B (, p. 17]: H®(Z) — H(Z) = H®(Z) , (2) 


H®(Z) = lim H,(Z, R; 0) 
R-A-0O 
_ 18 fay E(Z’) fax Kwa) epee 4 a 
Raety cling eae te is 
and 


which may be expressed in a way similiar to H(Z). 
Further 


which is to be determined by equation (2). 
We introduce the condition of periodicity setting 


E(Z + L) = E(Z) exp (78). (5) 


@ will be supposed to be real except when otherwise stated. It follows from 
(4) that: 


i } <0 , , 3 d vAG f , 
Hoa SP yf az n(z) | Hed) exp i x(Z' — Z) 
SS a aS —oo 
+L/2 + 0° 
ee hah se paz Ba +n © de SOA) exp in(Z’+nL —Z) 
~ ant, J, ( | v AC a 
iB +Lj2 +00 jie 
L , , 
it LY Ss. (ind) | : 
—L/2 —oo 
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If necessary we may assure the existence of the inner integral by the intro- 

duction of small losses (Im $B <0) which, in this case, remove the poles from 

the axis of integration exactly in the sense required by the Sommerfeld radia- 

tion condition [2, p. 196]. Consequently we can apply Porsson’s Theorem 

(3, p. 69] with the result: 

sg 2nnx+0 
ie 


1 Ms Jil A) expt (Z* Z) 


n=-—oco ” 


2nn+6\2 
papel eee 


This is exactly the expansion of H into spatial harmonics. The kernel 


| g"(Z, 2’; 8) -> nae ear pil 7 (6) 
has the following properties: 
(a) it is hermetic: AW ey AB 1 | ead WA Ag) ae 
(b) it is periodic in 6: g(Z, Z’; 0+ 22) = g(Z, Z'; 8) , (7) 
(c) g(Z, 2’; —8) = g*(Z, Z’; 8). 


The same properties will be assumed to hold for the kernel resulting from the 
field outside the electron beam. For example it may be shown by Fourier 
analysis that 

(a) for a ring-loaded waveguide (Figure 2) 


Figure 2 


o(® 4 Zz 0 os 1 Ji(Mm A) Yo rn B) pst J (Ym B) Yi (%A 
aes) Vy L Jo Vp A) Yo(% B) _ Talng B) Yo(% A (8) 
x exp 1(Z’ — Z) : Bin | 


where 


, 
aa 
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j 
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(b) for a disc-loaded waveguide (Figure 3) 


Figure 3 


Disc-loaded waveguide. 


a dn Jaltn-A) Yolttn B) — Jolitn B) YilunA) 


o(4) Lie 
5 (Z, Z’; C= fim yy. Jon A ) Yeti, 2) = Solis, B) Yo(un,A) 


Z 1 VAs 1 
X COSH 7 (5 =i 3) COSH It (>: aE 3) 


where 
je p* — (Sy: 6.= Vv n= Oe0, = 2," 0. 


ing modified Bessel functions if the arguments become imaginary. 


_ Thereby the occuring cylinder functions have to be replaced by the correspond- 


A similar transformation is now applied to H“)(Z). Equation (5) leads to 


the following Fourier representation: 
+ L/2 


2nn+6 
= = i naliy E (Z") d aZ" , es pom es idl 
465 . 
This is introduced into equation (3) yielding 
+L/2 : : 
€ ‘) N , eI, =e 
H(Z Jim ey iz? S) exp i (Z" + Z'—Z) 7 
. sine,Z ' 
ee te exp (—7 2’) 
k=1 k 
+L/2 
2nn+0 


=i 2h tim / EZ wm az" 18S expid "3 2) 


a Sin 4.2% feant+o ; 
x | az! - Caen exp i 7 1) aa 
0 
To justify the last step we have to impose the condition Rei 9< 0. We c 
thus calculate the inner integral. 


, sing, Z’ .(2nn+ 8 ie : 
fu aes expi (= "7 = -1)Z ~ e2 — [(2na+ CN aoe te 


an 
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Hence we have the representation 


: +L/2 i 
H®(Z) = +f [ e(Z, Z'; 0) E(Z') aZ’ 
-L/2 
where 
5 a yr expi(Z’— Z) (2na+4+ O/L 
gL: - ar oe ez— [(2n a+ 6)/L — 1]? (10) : 
which again, by analytic continuation, is valid for real (and complex) values | 
of 0. 
Setting 


g(Z, 2’; 8) = g(Z, 2’; 6) — g9(Z, 2’; 8) 
and presuming that within an elementary cell, the boundary condition is 
EZ) 20 ea 
our integral equation (2) becomes: 


+D/2 


i [e(Z, 2’; 0) — g(Z, Z'; 0) E(Z’) dz’ =0. (11) 


—pD/2 


§ 3. An Equivalent Variational Principle ([4] p. 79) 


Multiply (11) by E*(Z), integrate and divide by the square of the line 
integral of the conjugate-complex electric field: 


+ D/2 +D/2 
ii i [g(Z, 2’; 0) — g(Z, Z’; 0)] E*(Z) E(Z’) dz dz’ 
F(6) = 3 eo #0. (12) 
TOR EIEN SE Vat} 
—D/2 —D/2 


We calculate the changes in F(6) with small variations of E(Z) keeping 
6 = const: 


py — [[lalZ, 2 0) — OZ, 2°; )) (BE*(Z) BZ’) + BA(Z) 6E(Z’)) az az’ 
| [E(Z) az\2 Btn ic Pad | 
_ ffieZ, 2 6) — eZ, 2°; 6) E*(Z) E(Z’) az az’ 
([E(2) az4 
 [[l9E*Z) E(Z’) + E*(Z) dE (21) az az’ 
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As the kernel are hermetic it follows: 


dF (0) = Real : 6E* (Z) aZ ) az | [e(Z, 2; 0) — g(Z, 2"; 0)] E(Z’) dz’ 


if oe 


- ae ioe Zor Bz, Z'; 0)) E(Z’) dz" 
Se eee 


x 2 Real [8E*(Z) dZ f E(Z’) dz’ 


and by (11) ; 
| oF (6) =0. | (13) 


We can thus find a good approximation for F(6) if we take as E(Z) a function 
which is close to the rigorous solution of (11).. Such an estimation may be 
; provided by a suitable electrostatic problem. Hence E(Z) will be supposed to be 
a real, known function. 


§ 4. Derivation of an Electromagnetic Differ ence-Equation 


By defining voltage and influenced current with [1, p. 18] 


‘it follows from (12) 


where the admittance function G(@) has been introduced as: 


= een ee ee 
+D/2 +D/2 
[eZ, Z’; 0) E*(Z) E(Z’) aZ dZ’ 
2aBpA —Ddj2 -v/2 s 
G(O) =F HoaeDE (15) 
[ [ B*(@) E@) az az’ 
—D/2 —v/2 


ee eee es 
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By the argument following equation (13) G(#) may be supposed to be a known 
function of 6. It is easily proved by (7) that this function has the properties: _ 


, G(0) = G*(6);- G0 + 22) = G(6) | : 
and (16) : 
G(+.0) = G(0) as E*(Z)= E(Z). | : 


Although 6 was initially assumed to be real we may expand the range of : 
validity to complex values by analytic continuation. Thus G(§) will be consid- 

ered as a simply periodic meromorph function with the simple period 22. 

Because of the properties (16) G(#) is real for the following values of 
O=x+1y: 


6 Su Sp Op oOo ey 1 joie Ody 2 cate 


The zeros 6 of G(0) are located along these lines. Let + 0, be the zeros of 
G(0) within a period strip defined by — a < x <a (see Figure 4) and let them 
be arranged in such a way that 


0 < Imag 6, < Imag 0, < Imag 6, <... 0S Reh <a 
then it follows from (16) that : 
G( + 0) =0; +06M%=0,+2n0; k,|n| =0,1,2.... (17) 


O=x+iy —plane 


Figure 4 
The zeros of G(@) within a period-strip. 
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- The zeros 6” will now be supposed to be of the first order except when 6), = 0, a. 
_ Caucuy’s method of decomposition into partial fractions [5, p. 118] gives then: 


hy SOE No 1 a | 1, Hee ] 
a) G(0) 6% o=6% eee 6 — (6,+ 2” 2) 6,+2n% 
+ © 1 +00 1 | 1 

+d Re cay > lor@s 2nn) ew | 
where 
ert rats 1 , 0G(6) 

ray Vea = = 18 

SG NC C10) Ae 00 at (18) 


_ are the residues corresponding to the zeros + 6". 
The inner series may be easily summed up yielding 


ef. 


1 os 1 co i Ox sin Or 
G(8) G0) + 2 G’ (6x) [cote 2 cos @ — cos a | ; 9) 


_ This representation is valid if 


| AG) 
kar dba 
C 


where C is a contour extending to infinity in such a way that no zero of G(0) 
becomes a point of C (see Figure 4). 
Let 


Imag 0-1 < Yn < Imag 6,, n>N>0 


O= Ags 7 ; by + Real (6,7 — 9,); & = 0, 1, 2,... 


then a suitable contour may be defined by 


oo wet ST Se ee ee ee ee re eee 


6= xXtna+iry , = Vi SVE In 
= X+4Vy > —NN — % SX SX%trA, 
=—x—natiy, Vy SYS In» 
C= K—-1Vy_, —NU—%SX Sx HNN. 
: Now if 
E |G +i y)[>M +0, —co<y<too, 
i |G(x +1 %)|>M+0, Sapp Oa eH AN > Os (20) 
| yn = (n°) , Se ise re NO, 


=a ee Te hee 


| ities 
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it follows 
+n + %tnn e 
des 2 dx 
0< | eee ge LL aor bare) lates 
%é -Yn ~%y— nm 
7 es es cee ee d 
oT Gee aa ED ee ee 


We may check the conditions (20) in special cases but it seems to be extremely 
difficult to transfer them into general conditions for the respective kernel 
e(Z, 2’; 0). 

Equation (14) may finally be written 


I I I 6;, sin 9; ms 
iG(@) i G(0) rs pa 1 G’(8,) [cotg 2 cos 8 — cos 1 V. (21) 


If the electron beam is absent we have J = 0 and we get non-trivial solutions 
V + 0 for infinite sets of eigen-values 


6=0M=46,4+2n2. 
We shall fix our attention to one of these sets, say 
k=0: 0M =4+6,4+2n2. 


Treating the influenced current as a small perturbation we have 9 ~ 6") and 
thus very nearly 


ae iV. iis: 49 a": sin 6, 
iG(0) ~ iGO) * 766, [cotg 2 cosd— a | 
(22) 


sin O;, 


7 I 0, 
: = i G’(O,) [cots 2) eas Orme eek | 
To compare this with equation (1) we introduce 


V(n) = Vexp (=I), I(n) =I exp (—nT) 


with the result 


“= (2 eh 


= ee erie ot pe (23) 


5 
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Consequently we get: 


DE 0. 
anges 
G’(B0) ’ 


£ 


Rs i 0 7 6 sin 0, 
W=coa t+ GH cotg( 5") ia 2 G’(6,) [cote() re = me 


Hence we have proved the validity of equation (1) under quite general condi- 
tions and found a way to calculate explicitely the important parameters 
I’, Z, W. An example of this calculation will be given in a following paper. 


Acknowledgement 


The author wishes to express his appreciation to Prof. Dr. F. TANK for his 
‘encouragement and to Prof. Dr. E. StreEFEL for his helpful criticism. 


8 


A= wa/uy 


no 


= w b/Uuy 
b= 3-108 m sec * 


List of Symbols 


beam radius 

normalized beam radius 
waveguide radius 
normalized waveguide radius 
velocity of light in free space 


d slot width 

D=adluy normalized slot width 

e/m = —1,76- 1011 m? V—1 sec~? specific electronic charge 

E7(Z, R) axial component of the electric field intensity 


fr = (Px/ A) [1 + (be/A)*) 
‘Ho(Z, R; e) 


l 

L= wl/uy 
Pr 

Y 
R=or|[uUy 
Ww 
az 

‘A 

Z 

B= u/c 
1h 

T 

e€= Q/a 
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amplitudes of partial plasma waves 

azimuthal component of the magnetic field intensity 
length of elementary cell = slot spacing 
normalized slot spacing 

zeros of the Bessel function J,(*) 

radial coordinate 

normalized radial coordinate 

disturbing impedance 

axial coordinate 

normalized axial coordinate 

coupling impedance 

normalized electron velocity 

transfer constant for cold waves 

transfer constant for influenced waves 
normalized plasma frequency of infinite beam 


114 


on = olV/.+ lA)? 


& Ge 1078/36 2 A sec Vm 


fx 1202 VA-} 
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Myo = 42%-107~7 V sec A~* m- 


p= Bp? — x? 
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Q= V (e/m) Qo/€o 
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1 permittivity of free space 


wave impedance of free space 
phase shift per elementary cell 
axial separation parameter 


1 permeability of free space 


normalized plasma frequencies of finite beam 

| 
radial separation parameter 
space-charge density | 
angular frequency 
plasma frequency of infinite beam : 
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Zusammenfassung 


Auf Grund der Sechspoltheorie der elektrischen Nachrichtentechnik wurde an 
anderer Stelle [1, S. 23] eine elektromagnetische Differenzengleichung zweiter 
Ordnung abgeleitet. Dieses auf formalem Wege gewonnene Resultat ist fiir die 
Theorie der Traveling-Wave-Tubes von Bedeutung. Mit Hilfe der Maxwell- 
Lorentzschen Theorie wird nun gezeigt, dass die genannte Differenzengleichung 
unter viel allgemeineren Bedingungen zu Recht besteht. Der gleiche Weg fiihrt 
zu einer Methode, die auftretenden Parameter explizite zu berechnen. 
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Gibt es eine statistische Mechanik der Sternsysteme? 


Von Rupo.tF KurRtTH, Bern?) 


§ 1 


In der Literatur iiber Probleme der Stellardynamik finden sich vielfache 
‘Berufungen auf die statistische Mechanik. Fasst man sie alle zusammen, so 
scheint etwa die folgende Vorstellung von der Entwicklung eines Sternsystems 
verbreitet zu sein: 
« Das System sei als mechanisch abgeschlossene Gesamtheit von sehr vielen, 
-etwa  gravitierenden Massenpunkten gedacht. Befindet es sich im statistischen 
Gleichgewicht, so bleibt es darin, abgesehen von zufalligen kleinen Schwan- 
kungen. Befindet es sich nicht darin, so nimmt dem zweiten Hauptsatze der 
Thermodynamik entsprechend seine Entropie im Laufe der Zeit mit grosster 
Se brecheinlichkeit zu, es strebt das statistische Gleichgewicht an und erreicht 
es nach einer Zeitdauer von der Gréssenordnung der Relaxationszeit. Es 
besitzt dann eine Boltzmannsche Energie-, insbesondere also eine Maxwellsche 
Geschwindigkeitsverteilung. Ferner stellt sich Gleichverteilung der Energie ein. 
Ein kugelsymmetrisches System sollte dann wie eine isotherme Gaskugel auf- 
-gebaut sein, und zwar mit iiberall unendlich kleiner Dichte; denn andernfalls 
" miisste es eine unendliche Gesamtmasse besitzen. Den Prozess der Ausdehnung 
"kann man sich konkret als ein Entweichen der Sterne hoher Geschwindigkeit aus 
dem noch in Entwicklung befindlichen System vorstellen. Auflésungsgeschwin- 
digkeit und -zeit lassen sich an Hand der Maxwellschen Verteilung berechnen. 
Rotierende Systeme nehmen schliesslich eine konstante Winkelgeschwin- 
igkeit an. Im Milchstrassensystem ist weder dieser Zustand starrer Rotation 
och der der Gleichgewichtsverteilung erreicht; hieraus kann man sein Alter 


abschatzen. » 
Die Richtigkeit dieser soeben skizzierten Meinung soll im folgenden ge- 


_ priift werden. 
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§ 2 


Statistik kann man nur treiben, wenn sich ein dem Problem angemessener 
- Wahrscheinlichkeitsbegriff erklaren lasst. Fir das betrachtete mechanische 
- System von 7 Massenpunkten ist die Wahrscheinlichkeit als eine Mengenfunk- 
tion anzusetzen, die fiir die (messbaren) Teilmengen des (6 2)-dimensionalen 
Phasenraumes J” definiert ist und gewissen einschrankenden Bedingungen. 


ee 


_ = 


1) Astronomisches Institut der Universitat Bern. 
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unterliegt [8]1). Mit der Existenz einer solchen Wahrscheinlichkeit ist (im 

wesentlichen) die einer Wahrscheinlichkeitsdichte von folgenden Eigenschaften | 

logisch gleichwertig : 

1. Die Wahrscheinlichkeitsdichte ist eine im Phasenraum J” erklaérte Punkt- 
funktion. : 

2. Sie ist nichtnegativ. 

3. Sie ist normiert, das heisst, ihre Integration iiber den ganzen Phasenraum 
ergibt den Wert eins. 

4, Sie ist ein (im allgemeinen zeitunabhangiges) Integral der Bewegungsglei- 
chungen. 

Der letzte Punkt ist im wesentlichen eine unmittelbare Folge aus dem 
Liouvilleschen Satze. Ein System muss namlich fiir alle Punkte seiner Bahn- 
kurve im Phasenraum J” dieselbe (infinitesimale) Wahrscheinlichkeit besitzen — 
andernfalls wiirde der mathematische Begriff « Wahrscheinlichkeit » allen Bezug 
auf den allgemeinen Begriff verlieren. Nun bleibt nach jenem Satze der Inhalt 
eines Volumelements langs der Bahnkurven konstant. Also muss dasselbe auch 
fiir den Quotienten als Wahrscheinlichkeit und Inhalt gelten, das heisst fiir die 
Wahrscheinlichkeitsdichte. 

Wir k6énnen also sagen: Eine statistische Mechanik fiir eine bestimmte, 
durch die Bewegungsgleichungen gekennzeichnete Art von mechanischen Sy- 
stemen ist dann und nur dann méglich,wenn sich ein nichtnegatives normier- 
bares erstes Integral dieser Bewegungsgleichungen angeben lasst. Dieses Integral 
kann dann als Wahrscheinlichkeitsdichte gewahlt werden. 

Hierbei haben wir den Fall betrachtet, dass fiir den ganzen Phasenraum 
eine Wahrscheinlichkeitsverteilung definiert werden soll. Statt dessen kann 
man sich auch darauf beschranken, Wahrscheinlichkeiten tiberhaupt nur auf 
einer in J" gelegenen «invarianten Hyperflache» X' zu erklaren; es soll also Y’ 
dadurch definiert sein, dass ein oder einige erste Integrale der Bewegungsglei- 
chungen auf 2’ konstant sind. 

Inhaltlich liefert dieser neue Wahrscheinlichkeitsansatz nichts prinzipiell 
Neues. Auch fiir ihn bleibt die obige Bedingung im wesentlichen giiltig. In den 
praktisch wichtigen Fallen reduziert sie diese auf die Bedingung, dass »’ einen 
positiven endlichen Flacheninhalt besitzt. 


§ 3 


Im Anschluss an MAXWELL entwarf HECKMANN [5] eine statistische Mecha- 
nik gravitierender Massenpunktsysteme folgendermassen: Auf einer Hyper- 
flache 2’, die durch die Konstanz der bekannten elementaren Integrale erklart 
ist, wahlte er als (Flachen-) Wahrscheinlichkeitsdichte die entsprechende Ver- 
allgemeinerung der bekannten mikrokanonischen Verteilung. Dass die Nor- 


) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 125. 
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- mierungsbedingung erfiillt ist, setzte der Autor stillschweigend voraus. Das 
jedoch trifft wegen der besonderen Form des Newtonschen Potentials sicherlich 
nicht zu — wie wenigstens fiir den Fall » = 3 bereits PoINCARE gezeigt hat [9]. 
: Hier tritt also die grundlegende Schwierigkeit und, wie es scheint, Unmog- 
lichkeit einer statistischen Mechanik der Sternsysteme zu Tage: Es gelingt 
_ nicht, ein nichtnegatives normierbares erstes Integral der Bewegungsgleichun- 
- gen anzugeben. Das Gravitationspotent’al steigt bei wachsender Entfernung 
- zu langsam an, die Gravitationskrafte sind zu « weich», so dass die Flachen 2’ 
sich ins Unendliche erstrecken und unendlichen Inhalt besitzen. 


§ 4 

CHANDRASEKHAR und von NEUMANN versuchten den folgenden Entwurtf 
einer statistischen Mechanik der Sternsysteme. 
Sie betrachteten » Massenpunkte innerhalb einer Kugel K des realen 
_ Raumes #%, die den Radius R besitzt. Die Ortsvektoren der Massenpunkte 
 seien r,, ihre Impulse seien p; (¢ = 1, 2, ..., 2). Alle Massen werden als gleich 
angenommen. Fiir die Wahrscheinlichkeitsdichte 1,(r;, p,) des 7-ten Massen- 
_ punktes setzen die Autoren willkiirlich 

| 1 _f'_¢-? wenn 1,€ K, 
= 4 R3/3 mel (4, 1) 


T (Cis Pa) 
lo, wenn t;€ Rs — K, 


4 j = const > 0; und ferner fiir die Wahrscheinlichkeitsdichte T in 1”: 
T= [] vitti. Pi). (4. 2) 
i=1 


Hiermit bestimmen sie die Verteilungsfunktion und « mittlere Lebensdauer » der 
Gravitationskraft. An geeigneter Stelle wird dabei der Grenziibergang 


: nH 
n->oo, R->oo mit 4 R373 ==,constc>)0 


vollzogen. 
Jedoch: T ist fiir kein endliches n ein Integral der Bewegungsgleichungen, 


das System andert von Augenblick zu Augenblick seine (infinitesimale) Wahr- 


scheinlichkeit. lim T ist zwar ein Integral, ist namlich Null, ist aber gewiss 
Roo 


nicht normierbar und daher als Wahrscheinlichkeitsdichte ebenfalls unbrauch- 
bar. 


§ 5 
Wenn auch der soeben erdrterte Wahrscheinlichkeitsansatz nicht als zu- 
s doch erlaubt, ihn in ganz einfacher Weise aus- 


lassig erscheinen kann, so sel € 
h von einer andern Seite zu beleuchten. 


zuwerten und das Problem so noc 
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Die Autoren fragten nach der Verteilungsfunktion der Gravitationskraft 
(je Einheitsmasse), die die 7 Massenpunkte in irgendeinem festen Raumpunkt 
Y) erzeugen. Man kann nun eine « Statistik des Gravitationsfeldes in J”» ohne 
jede Rechnung aufstellen, wenn man statt der Kraft das Potential 


. Gm, 

Ce on 

auf die Einheitsmasse untersucht. (G bedeute die Gravitationskonstante, mM; 
die Masse des 7-ten Partikels.) 

Fir die Wahrscheinlichkeitsdichte im Phasenraum J” nehmen wir wieder den 
Produktausdruck (4. 2) als giiltig an, lassen aber den speziellen Ansatz (4. 1) 
fallen, indem 7,(t;, p;) nun irgendeine « verniinftige » Wahrscheinlichkeitsdichte 
bedeuten soll. Die Einzelpotentiale V; kénnen dann als unabhangige statistische 
Variablen angesehen werden. Nach dem Grenzwertsatz der Wahrscheinlichkeits- 
rechnung besitzt (unter geeigneten Voraussetzungen iiber die T;) das Gesamt- 
potential V eine annahernd GauBsche Verteilung: 

7 1 -(V-V)?/2 0 
pid Eason = (5. 2) 


wobei das Streuungsquadrat 


(V—V) a ee 


gesetzt ist. Querstriche bezeichnen Erwartungswerte. Diese beziehen sich zu- 
nachst auf die Verteilungsfunktion von V ; zum Beispiel: 


+co 
V=/Vo) dV, | 
hee canis (5. 3) 
ey [wv ~ V)2 w(V) av. 
Jedoch gilt offensichtlich fiir jede Funktion /(V): 
-- 00 
j= Hf iV) w(V) dV 


a LACH AMMEN PCCM AT Wen yh hater rer tty 
(5. 4) 


Durch die Schreibweise haben wir angedeutet, dass V in der oberen Zeile als 
unabhangige Veranderliche, in der unteren Zeile als Funktion der Vektoren 
Yo, YT, ..., Y, anzusehen ist. 
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Erklaren wir noch die Wahrscheinlichkeitsdichten #; und U durch 


u,(¥;) =f oiles, p,) 4p;, 
(pz) 
oo (5: 5) 


PU(ty,---.t) =/ seo f Teta yes tar Pav oo Pad dps, ..-, Bn = J Jule), 
(Hi) (Pn) ps 


so erhalten wir aus (5. 4): 
j= [-[ HV 0 tt] Wt ack. FRET REAL A; 
(t1) (tn) 


- ferner: 


V=D'V; 
i=l 
n iene 
= [r, — ti u,; (t,) dt;, 


3 


= 30, —0)?+ 3 Via Vd) Vi Vi 


(5. 6) 
| 


i=1 i+} 
=); (V;—Vi)* 
t=1 
“ —Gm 2 
=F fei — Maleate ae 
1=1 i 0 
(tz) 


das heisst: Die relative Streuung der Werte V um den Erwartungswert V (to) 
ist O(1 Vn), wird also bei hinreichend grosser Partikelzahl beliebig klein; das 
Potential V ist fiir die verschiedenen Konfigurationen (t,,T2,-++)T) nahezu 
konstant und nahezu gleich dem Erwartungswerte V (to). 

Und zweitens: n m ist die Gesamtmasse des Systems. Das Produkt 


o(r) =m u(t) 


i 
| 
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kann als seine Massendichte im realen Raum R* gedeutet werden — das heisst: 
Ist B irgendein Bereich des 3, so ist | 


f o(r) dr 


: 
der Erwartungswert fiir die in B enthaltene Masse. Fiir diese Massendichte 
gilt nun nach (5. 6): 


Vie) = —6 [era 
Di 


der Erwartungswert des Potentials ist also so zu berechnen, als ob die Masse des 
Systems iiber den realen Raum mit der Dichte g(r) stetig ausgebreitet ware. 

Eine geringe Modifikation des Gedankenganges ergibt ein ahnliches Resultat 
fiir das Potential, das je » — 1 Massenpunkte auf den jeweils -ten ausiiben. 
Das heisst also: Das « geglattete» Potential des ganzen Systems sollte tiber das 
von irgendeinem Einzelstern herriihrende Potential immer stark tiberwiegen. 
Insbesondere diirfte es keine « Doppelsterne» geben. Die Theorie miindet in die 
gewohnliche Stellardynamik ein. 

In den wirklichen Sternsystemen gibt es nun aber Doppelsterne, und zwar 
relativ recht haufig. Das wird im Rahmen jenes Ansatzes dann so gedeutet, dass 
wegen des geringen Alters zum Beispiel des Milchstrassensystems die Doppel- 
sterne noch nicht geniigend Zeit gehabt hatten, sich aufzulésen [2]. 

Uberhaupt erscheint die Existenz von Doppelsternen mit der Unméglich- 
keit einer statistischen Mechanik aufs innigste verkniipft: 

Doppelsterne scheint es (bei irgendeinem Gesetz der wechselseitigen Anzie- 
hung) dann und nur dann geben zu kénnen, wenn das Anziehungspotential 
zwischen je zwei Massenpunkten eine nach unten unbeschrankte Funktion ihres 

_Abstandes ist. Dann wird namlich bei geniigender Annaherung zweier Sterne 
das gegenseitige Potential absolut viel grosser als das Potential, das die tibrigen 
Partikel auf sie ausiiben. 

Sind umgekehrt die wechselseitigen Potentiale nach unten beschrankt, so 
k6nnen sie als nichtnegativ mit dem Minimum Null angenommen werden. Bei 
hinreichend grosser Partikelzahl n tiberwiegt dann aber das geglattete Gesamt- 
potential das gegenseitige Potential der beiden Massenpunkte so stark, dass 
diese im Laufe der Zeit wohl sicherlich getrennt werden. 

In diesem Falle lasst sich mittels des Energieintegrals gewiss eine Wahr- 
scheinlichkeitsdichte im Phasenraum J” konstruieren. 

Sind die Potentiale nach unten nicht beschrankt, so erstrecken die Flichen 
konstanter Energie sich bis ins Unendliche, und die Definition einer Wahr- 
scheinlichkeitsdichte wird héchstens in Ausnahmefillen méglich sein. — 


ee Ge 


enV 
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Zu unserer Untersuchung des Wahrscheinlichkeitsansatzes (4.2) sei noch 
angemerkt: Es ist nicht ohne weiteres méglich, in ihr das Potential durch eine 
Kraftkomponente zu ersetzen. Der Integralausdruck fiir das Streuungsquadrat 
wiirde namlich divergieren, wie man aus (5. 6) ohne Miihe abliest, indem man 
den entsprechenden Ausdruck fiir den Potentialgradienten bildet. 


§ 6 


Wir wissen, eine Wahrscheinlichkeitsvertei ung lasst sich im Phasenraume 
T sicherlich nicht definieren. Nehmen wir an, es ware doch méglich. Was 
wiirde daraus folgen ? 

Der Boltzmannsche Verteilungssatz kénnte dennoch nicht gelten. Denn fiir | 
ihn ist die Voraussetzung wesentlich, dass die Wechselwirkungsenergie gegen 
die iibrige mechanische Energie vernachlassigbar klein ist. Das trifft fiir Stern- 
systeme nun gerade nicht zu: Ihre gesamte potentielle Energie ist Wechsel- 
wirkungsenergie. 

Damit wird auch die Konstruktion der isothermen Gaskugel als eines Mo- 
delles fiir die Sternsysteme hinfallig — sie beruht wesentlich auf dem Boltz- 
mannschen Satze. 

Dagegen bliebe die Maxwellsche Geschwindigkeitsverteilung giiltig. Zwar 
ist sie eine spezielle Boltzmannverteilung, jedoch bleibt sie von dem obigen 
Einwand unberiihrt, weil es fiir die kinetische Energie keine Wechselwirkungs- 
glieder gibt. 

Auch der Gleichverteilungssatz wiirde fiir die kinetische Energie daher be- 
stehenbleiben. 

Dass er im Milchstrassensystem nicht erfiillt zu sein scheint, darf nach dem 
Friiheren aber nicht als Grundlage fiir irgendwelche Schliisse tiber das « Alter » 
das Systems benutzt werden. Im iibrigen hatten wir bei solcher Anwendung 
iiberhaupt erst noch zu klaren, in welcher Weise statistische Aussagen am kon- 
kreten System interpretiert werden kénnen [7] (8). 

Die Rechnungen iiber die « Auflésung» oder « Verdampfung» etwa eines 
Kugelsternhaufens [1] wiirden durchaus bedenklich erscheinen: Das Max- 
wellsche Gesetz der Geschwindigkeitsverteilung gilt nur approximativ [7] [8]. 
Bei jenen Rechnungen wird es nun aber gerade in Geschwindigkeitsbereichen 
benutzt, in denen eine geniigende Annaherung nicht verbiirgt werden kann. 

Wie sich tibrigens aus einem fundamentalen Satz der allgemeinen Mechanik 
von E. Horr [6] folgern lasst, kann die Auflosung eines Sternsystems «von 
selbst », das heisst ohne die Mitwirkung dusserer Krafte, nur ein ganz seltener 
Ausnahmefall sein — und nicht die Regel, wie man nach jenen Rechnungen an- 
nehmen miisste. (Damit der Hopfsche Satz anwendbar wird, muss durch pas- 
sende Abanderung des Gravitationsgesetzes fiir sehr kleine Entfernungen 7 die 
Singularitat des Potentials fiir 7 = 0 beseitigt werden.) 


122 Rupo.tr Kurta ZAMP 


Schliesslich: Die Berufung auf den zweiten Hauptsatz ware in der in §1 : 
wiedergegebenen Weise nicht statthaft. Jedoch kénnte das Gibbsche 1: -Theo- | 
rem angewandt werden. Da die genauere Erérterung dieses Punktes einen allzu 
breiten Raum einnehmen wiirde, verzichten wir hier auf sie und verweisen — 
auf [8]. | 


§7 


Schliesslich werde der Begriff der Relaxationszeit T noch etwas eingehender 
behandelt. Er wird zum Beispiel so definiert: Wahrend eines Zeitintervalls von 
der Lange T erreicht der mittlere Umsatz eines Sterns an kinetischer Energie _ 
deren Mittelwert — oder ahnlich. Gemeinsam ist allen Erklarungen: Gegen — 
Ende dieses Intervalls haben die Begegnungen eines Sterns mit anderen Einzel-_ 
sternen im Mittel einen merklichen Einfluss auf seine Bewegung gewonnen, so 
dass nun nicht mehr damit gerechnet werden kann, dass das in der Stellardyna- 
mik iibliche Modell einer im sechsdimensionalen Phasenraum yu strémenden 
Flissigkeit brauchbar bleibt. 

Es wird behauptet bzw. als selbstverstandlich angenommen: Innerhalb 
eines Zeitintervalls von der Lange T stellt sich das « statistische Gleichgewicht » 
her, sofern es vorher nicht bestand. Mit statistischem Gleichgewicht ist eine 
Boltzmannsche bzw. eine Maxwellsche Verteilung gemeint. 

Unsere vorangehenden kritischen Erérterungen haben jedoch wohl hin- 
reichend dargetan: Diese Behauptung k6nnte beziiglich der Maxwellschen 
Geschwindigkeitsverteilung nur dann als begriindet, ja tiberhaupt als sinnvoll 
angesehen werden, wenn im Phasenraum J" sich eine verniinftige Wahrschein- 
lichkeitsdichte angeben lasst. Und das allgemeine Boltzmannsche Verteilungs- 
gesetz wiirde auf keinen Fall gelten. 

Die iibliche Berechnungsart dieser Relaxationszeit benutzt die (sich ja 
eigentlich ausschliessenden) Grundvorstellungen von Partikel und Kontinuum 
zugleich, setzt gleichférmige Dichte und Maxwellsche Geschwindigkeitsvertei- 
lung voraus und beriicksichtigt nur die Begegnungen je zweier Sterne — im 
ganzen ein recht zweifelhafter Ansatz (abgesehen von dem obigen Kardinal- 
einwand), CHANDRASEKHARS Rechnung fiir kugelsymmetrische Systeme [1] 
zeitigt das paradoxe Ergebnis: Je mehr Sterne innerhalb einer Kugel von 
festem Radius, um so grésser T — um so kleiner also der Einfluss der Begeg- 
nungen. Das entspricht ganz den Folgerungen, die in §5 aus der Statistik 
(4. 1), (4. 2) gezogen wurden. 

Mittels jener Statistik versuchten CHANDRASEKHAR und von N EUMANN, den 
Mehrfachbegegnungen Rechnung zu tragen. Den grundsitzlichen Bedenken 
von §4 und §5 fiigen wir einen. empirischen Hinweis von FRICKE hinzu [4]. 
Demnach liefern entsprechende und entsprechend hergeleitete Formeln fiir die 
Relaxationszeiten in Elektronenplasmen Werte, die um Faktoren bis zu 103 
gegentiber den gemessenen zu gross sind. 
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§ 8 


_ Als Kern der Erklarung in § 7 fiir die Relaxationszeit T heben wir heraus: 
_ Gegen Ende eines Zeitintervalls von der Lange T haben die wechselseitigen 
_ Begegnungen der Sterne auf ihre Bewegungen einen solchen Einfluss gewonnen, 
_ dass das stellardynamische Kontinuummodell mit seinem geglatteten Potential 
_fernerhin nicht mehr als eine zureichende Naherungskonstruktion angesehen 
_ werden kann. 

Wir geben eine einfache Abschatzung von T. Sie ist zwar hypothetisch, er- 

scheint aber doch wohl als hinlanglich begriindet: 

Das Kontinuummodell wird durch die folgenden Grundgleichungen be- 
_ schrieben: 


of , yr Of OE bo ea: 
ae, ise sige Hs 5) ” 
+00 
0 hea: f(u, x, t) du, dug dus, (8. 1) 
3 
02V 
py Ox, =42Go 
j=l 
Hierbei bedeuten : 
sy die Zeit, 
x5 (i = 1, 2, 3) rechtwinklige kartesische Koordinaten im gewohnlichen 
Raum 3, 


u die zugeordneten linearen Geschwindigkeitskomponenten, 
& den Vektor (1, %2, %3), 

u den Vektor (#1, U2, 3), 

o(x,t) die Massendichte des Systems im ‘R®, 

}(u, x, t) seine Massendichte im sechsdimensionalen (uv, x)-Raum, ' 

V(x, t) das vom System erzeugte Gravitationspotential auf die Einheitsmasse. 
Wir fiihren nun dimensionslose: Variable ein; sie seien durch ~ gekenn- 


2. 


zeichnet : aan, work, 
pe GK, u,; = «* VG 0* “,, t= Vee® t 
(Ae es CMY 
x3 G VG oF 


x«* und e* bedeuten dimensionsbehaftete Konstante; x* eine fiir das jewel- 
lige Modell charakteristische Lange (Durchmesser, Radius oder dergleichen), o* 
eine charakteristische Dichte, etwa die mittlere. 
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Fiihrt man nun in die Grundgleichungen (8. 1) die dimensionslosen Veran-_ 
derlichen ein, so ergeben sich fiir diese ganz entsprechende Gleichungen mit : 
Getz jeder dimensionslosen Lésung ih gehért also eine zweiparametrige — 
Schar dimensionsbehafteter Lésungen /. : 

Der dimensionslosen Zeitdauer 1 entspricht die « Zeitskala » ¢* des Systems: — 


L pa 1,5 -107 
VG @* Vor 
wenn die Dichte in Sonnenmassen je Kubikparsek gemessen wird. 


Wenn wir nun sagen, das zeitliche Verhalten eines Sternsystems lasse sich | 
tiberhaupt durch das Kontinuummodell beschreiben, so kann das doch nur 


Hel 


[ Jahre], 


bedeuten: f bzw. f ist ee brauchbare Anndherung fiir ein Zeitintervall minde- 
stens von der dimensionslosen Gréssenordnung 1. Denn andernfalls wiirde ja 
die Anwendbarkeit des Modells unmittelbar nach dem Zeitnullpunkt aufhéren. 

Das heisst: Ist fiir ein Sternsystem das Kontinuummodell tiberhaupt an- 
wendbar, so ist die « Relaxationszeit» mindestens von der Gréssenordnung 
seiner charakteristischen Zeitskala. In Zeichen: 


1 10? 


Das ist exakt, weil ja fast eine Tautologie. Es bleiben jedoch noch zwei 
Punkte zu klaren: 

Erstens: Unter welchen Bedingungen ist denn die Anwendung jenes Modells 
zulassig? Uns erscheint die Vermutung plausibel: fiir alle Anfangszustande 
des Massenpunktsystems, die sich mit hinreichender Annaherung durch eine 
kontinuierliche Anfangs-« Verteilung» /(u, x, 0) darstellen lassen. 

Und zweitens: Wir haben fiir den Parameter die mittlere (Anfangs-) Dichte 
gewahlt. Man kénnte doch statt dessen auch zum Beispiel die maximale An- 
fangsdichte oder das 10?°-fache oder auch das 10-2°-fache dieses Wertes neh- 
men. Wie ist die Auszeichnung der mittleren Dichte zu begriinden ? 

Man fithre in die Bewegungsgleichungen, die zur ersten Gleichung (8. 1) 
gehoren: 


Ne 
heats 
ange o(y, t) 
ae vee we ly—«| dy, dy, dys 


J 3 


die dimensionslosen Veranderlichen ein. Die dimensionslosen Beschleunigungen 
werden dann (dem Betrage nach) durchschnittlich in der Gréssenordnung 1 
legen, wenn o* die mittlere Dichte bedeutet — ausser in ganz kiinstlichen Aus- 
nahmefallen. Das rechtfertigt die besondere Wahl fiir o*. — Die theoretische Er- 
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fahrung an mechanischen Modellen bestatigt die Richtigkeit dieser Uberlegung 
wohl ausnahmslos. 


Beschluss 


Auf die in der Uberschrift gestellte Frage haben wir mit einem klaren Nein 


- za antworten: Bis jetzt existiert eine statistische Mechanik der Sternsysteme 


nicht, und es ist auch nicht abzusehen, wie sie méglich sein k6énnte. 

Der Grund hierfiir liegt im Charakter der Gravitation: Sie ist nicht stark 
genug, das System auf jeden Fall zusammenzuhalten. Eine statistische Mechanik 
der Gase ist nur méglich, weil das Gas zum Beispiel in einem Gefass einge- 
schlossen ist oder durch die Schwerkraft der Erde am Entweichen verhindert 
wird. Man kann sagen: Fiir die statistische Mechanik ist das Gefass ein wesent- 
licher Teil des Gases, und dieses Gefass fehlt bei Sternsystemen. 

Es ist dies im Grunde dieselbe Schwierigkeit, die uns in der Himmelsmecha- 
nik bei der Aufstellung von Stabilitatskriterien begegnet: Die Flachen kon- 
stanter Energie und alle andern aus ihnen abgeleiteten « Hillschen Grenz- 
flachen » erstrecken sich ins Unendliche. 
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(klassischen) statistischen Mechanik. (Noch 


de la mécanique céleste, Bd. 3 (Gauthier- 


Summary 


It is shown that, on account of the ““softness”’ of the force of gravitation, 
statistical mechanics cannot be applied to star systems, and that therefore all 


results basing upon it have no fundament. 


(Eingegangen: 15. Marz 1954.) 
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Die Erregung reiner Eigenschwingungen von 
Flugzeugfliigeln’) 
Eine Anwendung der Theorie der Integralgleichungen 


Von ERNST SCHULTZE, Altenrhein?) 


1. Einleitung 


In der Flugzeugindustrie ist es tiblich, nach Fertigstellung eines neuen Flug- © 
zeugtyps vor der Flugerprobung den sogenannten Standschwingungsversuch 
durchzufiihren. Dabei werden mit elektrodynamischen Erregern sinusférmig 
variierende Krafte an einigen Stellen des Flugzeuges, insbesondere am Fliigel 
und am Leitwerk, angebracht. Der Zweck des Standschwingungsversuches ist 
die Ermittlung der Eigenfrequenzen und Eigenschwingungsformen von Fliigei 
und Leitwerk, deren Kenntnis dann fiir die Berechnung der kritischen Flatter- 
geschwindigkeit wichtig ist. 

Die Bestimmung der Eigenfrequenzen ist relativ einfach: Man variiert die 
Frequenz der erregenden Krafte so lange, bis eine Resonanz (Amplitudenmaxi- 
mum) eintritt. Schwieriger ist es, die zu dieser Eigenfrequenz gehérige Eigen- — 
schwingungsform (Amplitudenverteilung am Fliigel bzw. Leitwerk) zu bestim- 
men, besonders, wenn eine zweite Eigenfrequenz in der Nihe liegt. Denn obwohl 
das Flugzeug mit nur einer Eigenfrequenz erregt wird, wird die Amplitudenver- 
teilung ausser der zu dieser Eigenfrequenz gehérigen Eigenschwingungsform 
auch Anteile der andern Eigenschwingungsformen enthalten. 

Es entstehen deshalb die Forderungen 
1. nach einem Kriterium, ob in einem gegebenen Falle eine reine Eigenschwin- 

gungsform vorliegt oder nicht; 

- nach einer Methode zur Elimination méglichst vieler unerwiinschter Eigen- 
schwingungsformen, das heisst zur Herstellung einer méglichst reinen Eigen- 
schwingungsform. 

Es wird in dieser Arbeit gezeigt, dass 1. eine reine Eigenschwingungsform 
dann und nur dann vorliegt, wenn alle Punkte mit derselben Phase schwingen, 
und 2. um so mehr unerwiinschte Eigenschwingungen eliminiert werden kénnen, 


bo 


je grésser die Zahl der Erreger ist, wobei die Elimination durch Verandern der 


1) Diese Arbeit, von der bereits ein Auszug im Journal of the Aeronautical Sciences, August 
1954, erschienen ist, wurde im Anschluss an Untersuchungen geschrieben, die der Autor in den 
Flugzeugwerken Altenrhein gemacht hat. : 

*) Flug- und Fahrzeugwerke AG. 
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Kraftamplituden der Erreger vorgenommen werden kann und zum Teil auch 
durch Verschieben der Angriffspunkte der Erreger. 7 
Die Ergebnisse dieser Arbeit sind physikalisch sehr anschaulich und ver- 
mutlich gerade deshalb dem Praktiker zum Teil schon bekannt. Es herrschen 
aber sehr oft verschwommene Begriffe auf diesem Gebiet, was daher kommen 
mag, dass offenbar noch nirgends in der ohnehin sparlichen Literatur tiber den 


- Standschwingungsversuch die grundlegenden Zusammenhange klar dargestellt 


worden sind. Insbesondere ist die eigentlich naheliegende Anwendung der 
Theorie der linearen Integralgleichungen auf den Standschwingungsversuch 
offensichtlich neu. Es wird jedoch kaum méglich sein, ohne die Theorie der 
Integralgleichungen eine auch nur annahernd so einfache und iibersichtliche 
Herleitung der wesentlichen Gesetze zu geben. 

Der Flugzeugfliigel sowie die Leitwerkflachen konnen idealisiert werden als 
zweidimensionale ebene Flachen, die eine eigene Steifigkeit und eine Massen- 
verteilung besitzen und die langs eines Stiickes ihres Randes (am Rumpf) ein- 
gespannt sind. Sie fallen also unter den in der mathematischen Physik tiblichen 
Begriff «Platte». Eigenschwingungen von Platten sind zwar in der mathema- 
tischen Physik schon oft behandelt worden [1]1), jedoch nicht unter den ftir den 
Standschwingungsversuch geltenden Gesichtspunkten, welche oben erwahnt 
wurden. 

- Ausserdem lasst sich mit nur unwesentlichem Mehraufwand eine Art von 
Materialdampfung beriicksichtigen, die — im Gegensatz zu der sonst bekannten 
geschwindigkeitsproportionalen Dampfung — proportional zur Amplitude und 
in Phase mit der Geschwindigkeit ist. Eine solche Dampfung ist dann vorhan- 
den, wenn das Material im Spannungs-Dehnungs-Diagramm eine ellipsenfér- 
mige Hysteresiskurve besitzt [2] [3]. Die Beriicksichtigung dieser Material- 
dampfung ist besonders einfach, wenn man fiir die harmonische Zeitabhangig- 
keit die komplexe Schreibweise einfithrt ; man braucht dann namlich nur den 
Elastizitatsmodul des Materials mit der komplexen Zahl (1 + ¢ g) zu multiphi- 
zieren, wobei g eine Materialkonstante ist (fiir Metalle: g 0,025 hy0,08). 
Charakteristisch fiir diese Art von Dampfung ist, dass die Resonanzfrequenz 


von der Dampfung unabhangig ist. 
2. Bezeichnungen 


— kartesische Koordinaten in der Fliigelebene ; 


Hy, %. 
P(X,-%z} — Punkt auf dem Fliigel; ; 
y(P, t) — Verschiebung des Punktes P in vertikaler Richtung bei einer 


Deformation des Fliigels als Funktion der Zeit ¢; 
dP = dx, dx, = Flachenelement beim Punkte P; 
m/(P) = Massenverteilung ; 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 135. 
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m(P) dP = Masse des Flachenelementes dP; ; 
O(P, Q) = elastische Einflussfunktion (Greensche Funktion) 
= Deformation bei P, wenn Einheitskraft in vertikaler Richtung 
bei Q; 
v = Kreisfrequenz der Schwingung; 
y(P, t) = y(P) e'"; (1) 
y(P) = vertikale Amplitude des Punktes P, im allgemeinen eine kom- 


plexe Zahl, deren Betrag die Amplitude der Schwingung und 
deren Argument die Phase der Schwingung gegeniiber der 
Kosinusschwingung bedeutet ; 


g = Konstante der Materialdampfung; 

A = Eigenwertparameter; 

Q1, Qe, --., Qn = Punkte des Fliigels, an denen die  dusseren Erregerkrafte 
K,(t) angreifen; 

Ky (t), K,(é), ..., K,(¢) = vertikale Erregerkrafte, von denen vorausgesetzt wird, 


dass alle in Phase mit der Kosinusschwingung sind, also 


K,(t) = K, e*”* *(K;, reell, 7 =1, 2,..., n) , (2) 
ferner: 
p(P) = y(P) Vm(P) , (3) 
K(P, Q) = d(P, Q) Vm(P) m(Q) , (4) 
K; 
[paieeye Al 
7 ¥m(Q;) ' ©) 
AP) = (Lig) 3/8, (P,Q), 6) 
an = 37 Ky yu (Q) = 3 bs oe (Q) 7) 


3. Die Integralgleichungen und deren Lésungen 


Die infinitesimale d’Alembertsche Tragheitskraft dT in einem Punkte Q 
des Fliigels ergibt sich wegen (1) zu 


dT (Q, t) = »? y(Q, t) m(Q) dQ. 
Diese ergibt in einem Punkte P die Deformation 


dy(P, t) = d(P, Q) dT(Q, t) = »* 6(P, Q) m(Q) y(Q, t) dQ. 


er 


eT LO eT Lee 
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Die Integration iiber samtliche Punkte Q des Fliigels ergibt, nach Kiirzung des 
Faktors e'”', die folgende Integralgleichung fiir die Amplitudenverteilung y(P) 
des ohne dussere Krafte und ohne Dampfung schwingenden Fliigels: 


y(P) = 7 [ OP, Q) m(Q) ¥(Q) dQ. 8) 

Die Beriicksichtigung der Materialdampfung, wie sie in der Einleitung erwahnt 

wurde, erfordert die Multiplikation der Einflussfunktion mit der komplexen 
Zahl (1 — 7 g). Damit lautet dann die Integralgleichung 

y(P) = v8 (L = ig) [ 8(P, Q) m(Q) ¥(Q) ae. (9) 


Die ausseren Erregerkrafte liefern noch ein zusatzliches Glied zu dieser Inte- 
gralgleichung; sie erzeugen namlich im Punkte P die Deformationsamplitude 


> 6(P,Q;) K; im Falle ohne Dampfung, 
jot 

(1—ig) ¥” 6(P, Q;) K; im Falle mit Dampfung, 
j=1 


welche zu der von den Tragheitskraften erzeugten Deformationsamplitude 
addiert werden muss, um die totale Deformationsamplitude zu erhalten. Wenn 


- wir nur den Fall mit Dampfung ins Auge fassen, so lautet die (jetzt inhomogene) 


Integralgleichung : 
y(P) = 98 (1-1 9) [A(P, @) mQ) 1) €O-+ LF) SAPO) K-10) 


Alle drei Integralgleichungen (8), (9) und (10) haben den Kern 0(P, Q) m(Q), 
der nicht symmetrisch ist. Formuliert man sie aber fiir die in Gleichung (3) 


; definierte Funktion g(P), so lauten sie im Falle der freien Schwingung unter 


Beriicksichtigung der Bezeichnung (4): 
pP) = 2] K(P,Q) #10) 40, (11) 


~ wobel 


y? im Falle ohne Dampfung, | 


(12) 


| y2(1—ig) im Falle mit Dampfung | 
und im Falle der erzwungenen Schwingung mit Dampfung wegen (6) und (5): 


p(P) — 4/ KP, Q) (0) 40 = HP), (13) 


ZAMP. V1/9 
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wobei ee (14) 


Der Kern K(P, Q) dieser Gleichungen ist symmetrisch, weil die Einfluss- 
funktion 6(P, Q) symmetrisch ist (Satz von MAXwELL). Ausserdem ist er sogar 
positiv definit, was, kurz ausgedriickt, daraus folgt, dass man beim Aufbringen 
irgendeiner Belastung auf den vorher unbelasteten Fliigel immer eine positive 
Formanderungsarbeit leisten muss. 

Uber solche lineare Integralgleichungen zweiter Art existiert eine ganze 
Menge von Satzen [4]. So ist zum Beispiel die homogene Gleichung (11) nur 
lésbar fiir gewisse Werte von A, die sogenannten Eigenwerte. Diese sind wegen 
der Symmetrie des Kernes K(P, Q) reell, und weil er positiv definit ist, sogar 
positiv. Sie seien, der Grésse nach geordnet, bezeichnet mit 


MRO pee 


Zu jedem Eigenwert A; gehort als Losung der homogenen Gleichung eine reelle 
Eigenfunktion ,(P). (Den fir einen Flugzeugfliigel wenig wahrscheinlichen 
Fall mehrfacher Eigenwerte wollen wir hier tibergehen.) Es sei noch erwahnt, 
' dass die Eigenfunktionen zueinander orthogonal sind und dass man sie als 
normiert voraussetzen kann, also 
1 firi=k, 
J dP) gelP) AP = bie = : (15) 
0° fur? ef . 


Weil die Eigenwerte A, reell und positiv sind, folgt die anschauliche Tat- 

sache, dass 

a) im Falle der freien Schwingung ohne Dampfung die Gréssen », = Vax reell 
sind, die Schwingung also rein sinusférmig ist; 

b) im Falle der freien Schwingung mit Dampfung wegen (12) », nicht reell, 
sondern komplex ist: 


eS / 8 
Vi —_ | 1 a 2 == } Ax (1 a L | . (16) 


Fir den Bewegungsablauf hat man daher wegen (1): 


girkt = piVaxt eb VAnt/2 
das heisst, die Schwingung ist gedimpft, wie es ja sein muss. 

Weil ferner mit der Eigenfunktion g,(P) auch die Funktion yx(P) reell ist, 
folgt, dass bei einer Eigenschwingung der ganze Fliigel mit derselben Phase 
schwingt. 

Das Wesen des Standschwingungsversuches erfordert es, die Lésung der 
inhomogenen Gleichung (13) zu den Lésungen der homogenen Gleichung in 


ea es 


a eo ee 


Te CO ee a a 
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Beziehung zu bringen. Wir ziehen deshalb die bekannte Schmidtsche Reihe 
fiir die Lésung g(P) von (13) heran: 


oP) = 1(P) +A S75 oul), (17) 
wobei 
f=] HP) GP) aP. (18) 


v 


Weil 2 bei Anwesenheit von Dampfung nach (14) komplex ist, werden die 
Nenner A, — A nie Null, entsprechend der bekannten Tatsache, dass die Reso- 
nanzamplitude bei Materialdampfung endlich bleibt. Die Gleichung (17) lasst 
sich noch vereinfachen, wenn man /(P) und f/, durch die Eigenfunktionen und 


_ Eigenwerte ausdriickt. Dies wird dadurch erméglicht, dass man in der Formel 


(6) die nach dem Satz von MERCER gleichmassig konvergente Entwicklung 


- | @el(P) ox 
RProy=s y° Pe jel) 
k=1 K 


einsetzt. Die Integration in Gleichung (18) ist dann gliedweise gestattet, und 


die Formel (15) kann angewendet werden. Damit ergibt sich aus (17) schliess- 


lich die gleichmassig konvergente Reihe: 


o(P)= (1-58) Sy wulP) (09) 


wobei 


oder, indem man A, = 17 setzt und die Formeln (14) und (3) und (5) beriick- 
sichtigt : 


yPlm Se etresy HP) (20) 


2 ig pe 
_ + 1g % 
oa v*) + 18 ME 


wobei 


ay = i K; ¥e(Q3) « 
j=1 


Hierbei wird y,(P) durch die aus (15) folgende Bedingung 
/ [y(P) 2 m(P) dP = 1 (21) 


normiert und mit der Dimension L'”? K™ 1/2 T-1 behaftet. 
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Die Formel (20) sagt aus, dass bei Erregung des Fliigels mit sinusformigen 


Kraften gleicher Frequenz und gleicher Phase im allgemeinen sémtliche Eigen- | 
schwingungen (im eingeschwungenen Zustand) im Schwingungsbild des Fli- — 


gels vertreten sind, und zwar mit verschiedener Phase. Allerdings wird, wenn 


die Erregerfrequenz y gerade gleich einer Erregerfrequenz y,, ist, der Koeffi- { 
zient des zugehérigen y,,(P) wegen der Kleinheit von g bedeutend grésser als _ 
die andern Koeffizienten, so dass dann im allgemeinen die betreffende Eigen- — 
schwingung vorherrscht. Es kann jedoch unter Umstanden passieren, dass der _ 


Koeffizient a,, viel kleiner ist als die anderen Koeffizienten a,,, (k + m), so dass 


I 
| 


die Eigenschwingung y,,(P) doch nicht richtig zur Geltung kommt. Ausserdem | 


werden die Knotenlinien einer vorherrschenden Eigenschwingung durch die 
andern Eigenschwingungen verwischt. Es lohnt sich deshalb, die Méglichkeit 
der Elimination unerwiinschter Eigenschwingungen zu untersuchen. 


4. Erregung einer reinen Eigenschwingung 


Zunachst wollen wir beweisen, dass, wenn der erregte Fliigel in jedem 
Punkte mit derselben Phase schwingt, eine reine Eigenschwingung vorliegt. Aus 
Formel (19) und (14) folgt wegen A, = 07: 

Pls See ee at Pi 
p( ) a (v}— yp?) ae ig vi @x(P) (22) 
Das nach 'Voraussetzung von P unabhangige Argument der komplexen Zahl 
g(P) sei mit « bezeichnet. Dann gilt: 


Im [(P)] — tga Re[y(P)] = 0, 
also 
$7 ay [veg + (oR — v2) tea] wi 
z (ot — v8)8 + gt of g(P) =0. (23) 
Die Koeffizienten von @,(P) in dieser unendlichen Reihe sind von P unabhan- 
gig. Weil die Reihe (22) gleichmassig konvergent ist, ist es auch die Reihe (23). 
Infolgedessen darf man sie, auch nach Multiplikation mit der stetigen Funktion 
yi(P), gliedweise integrieren. Dies ergibt unter Beriicksichtigung von (15): 
a; [Mg + 0? —»*)tga]=0 fir7=1,2,3,.... 
Der zweite Faktor kann bei festem vy und « héchstens fiir ein 7 verschwinden 
(sagen wir fiir 7 = m), und zwar auch dann nur, wenn die Bedingung 


tga = — 2°", (24) 


j 


a 


¢ 
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: fiir die Phase « erfiillt ist. Somit hat man 


7 


@=0 tiri +m, (25) 


und die Reihe (20) degeneriert zu einem einzigen Gliede: 


: 4\(P) a aE 4m Yn F) hw ce 2 a an © es hae) ) 


A eae ig ye 
(Vin v ) i ts § Vm g Pm 


das heisst, der Fliigel schwingt in einer einzigen Eigenschwingung; was zu 
_ beweisen war. 


: 


Obwohl man theoretisch mit endlich vielen Erregern nicht unendlich viele 


_ Eigenschwingungen eliminieren kann, ist es doch in der Praxis méglich, reine 


Eigenschwingungen zu erhalten. Denn beim Standschwingungsversuch interes- 
siert man sich nur fiir die untersten paar Eigenschwingungen, so dass man bet 


_ der Erregerfrequenz nie iiber eine gewisse Schranke hinausgehen wird. Aus der 


Konvyergenz der Reihe (20) folgt dann, dass von einer gewissen, geniigend hohen 
Ordnung an keine Eigenschwingung mehr storen wird. 

Mit  Erregern kann man mindestens » — 1 Eigenschwingungen eliminieren, 
denn man hat die Bedingungen a, = 0 fiir k + m, welche ausgeschrieben fol- 


_ gendermassen lauten: 


ky Vi + Ky Vi2 re +.K, Vin == (0 
Ky Van Rove to + an Ue 


I r r = 
Ky Vm—1,1 os Ke Vm—1,2 = ai ap I, Vin ne 0 ) 
ky Vm+11 ie Ky Vm +1, 2 | ie a Ke Vir: 19g oe Os 


Ky Vni Ss Ky Vne2 Siang oe I Vnn = 0, 


ky Vn +1,1 aM Ky Vn+1, 2 he EA By, Yn+1,n 0, 


wobei y;,(Q;) = Yui gesetzt wurde. Setzt man zunichst einmal Q; als gegeben 
voraus, so hat man ein lineares homogenes Gleichungssystem fiir die Unbe- 
feannten A, io, «- Ky Nach den Regeln der Theorie der linearen homogenen 
Gleichungssysteme existiert eine bis auf einen konstanten Faktor bestimmte 
Lésung (Kj, ..., K,), welche mindestens 7 — 1 der obigen Gleichungen befrie- 
digt. Dies bedeutet, dass man bei festgehaltenen Erregerstellen Q; durch geeig- 
nete Wahl der Kraftamplituden Ky, ..., K,, mindestens 1 — 1 Eigenschwingun- 
gen eliminieren kann, 


ai 
» 
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Hat aber das Gleichungssystem der ersten 7 —1 Gleichungen einen Rang 
[= Rang der Matrix Y = (y,,)], der kleiner als 7 — 1 ist, so lassen sich weitere 
Gleichungen erfiillen. Dadurch werden dann mehr als » — 1 Eigenschwingungen — 
eliminiert. Eine Verkleinerung des Ranges lasst sich aber immer durch geeignete 
Wahl der Erregerstellen Q; erreichen, zum Beispiel dadurch, dass man samtliche | 
Q,; in die Knotenlinie einer Eigenschwingung legt; dann wird diese Eigen- 
schwingung von selbst eliminiert. Um diese Tatsache auszuniitzen, wird man | 
also die Erregerstellen in die Schnittpunkte von Knotenlinien verschiedener 
unerwiinschter Eigenschwingungen legen. 

Es leuchtet deshalb ein, dass man mit wenigen Erregern praktisch gentigend | 
reine Eigenschwingungen erzielen kann, wenn man sowohl die Grésse der 
Kraftamplituden als auch die Angriffspunkte der Erreger variiert. 

Abschliessend sei noch der Fall einer stetigen Verteilung der Erregerkrafte 
angefihrt, der dann von Interesse ist, wenn man den Aufwand sehr vieler Er- 
reger nicht scheut. Dabei werden wir auf eine andere Methode zur Herstellung 
reiner Eigenschwingungen gefiihrt, welche praktisch schon verwirklicht worden 
ist [5]. 

Die Erregerkraft pro Flacheneinheit sei bezeichnet mit K(Q), und entspre- 
chend zu Formel (5) sei 


k(Q) = d 
(Q) Vm(Q) (5a) 
Formel (6) geht dann iiber in 
((P) = (Lig) | K(P, Q) HQ) dQ (6a) 
und Formel (7) in 
ay = | BQ) pulQ) 40 = [ K(Q) (0) a. (7a) 


“ 


Deshalb lautet die Bedingung (25) fiir die Erregung einer reinen Eigenschwin- 


gung, dass k(Q) zu allen @,(Q) orthogonal sein muss ausser zu @,,(Q). Nimmt 
man daher 
k(Q) proportional zu @,,(Q) , (26) 


so ist a, = 0 fiir k + m. Aus (26) folgt wegen (5a) und (ore 
K(Q) proportional zu y,,(Q) m(Q) (27) 


als Bedingung dafiir, dass nur die m-te Eigenschwingung erregt wird [6]. 

Gehen wir wieder zu dem praktisch allein realisierbaren Fall endlich vieler 
(aber sehr vieler) Erreger itber, so wird naherungsweise dann nur die m-te 
Eigenschwingung erregt, wenn 


unabhangig von 7 (28) 


; 
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ist. Dazu sei erwahnt, dass Lewis und WRrisLey [5], die offenbar durch prak- 
tische Erfahrungen mehr gefiithlsmassig auf die Bedingung (28) gekommen sind, 
eine Apparatur entwickelt haben, die fiir jeden Erreger eine zur Verschiebungs- 
amplitude y(Q,;) proportionale Kraftamplitude K;, liefert. Es lassen sich damit 
die reinen Eigenschwingungen beinahe automatisch herstellen. Ausserdem 
-ertibrigt sich das Verschieben der Angriffspunkte Q; wegen der Vielzahl der 
_ Erreger. 
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Summary 


The steady-state response of airplane wings due to externally applied sinusoi- 
dal shaker forces is theoretically examined by means of the theory of linear 
integral equations of the second kind. Damping is accounted for, and emphasis 
is laid on the problem of exciting pure natural modes. The following items are 
proved: 
(1) Equal phase at every point of the vibrating wing is shown to be a criterion 
for the presence of a pure natural mode. 

(2) » shakers allow to eliminate at least n — 1 natural modes by variation of the 
exciting amplitudes. By displacing the shakers still more natural modes can 
be eliminated. 

The experimental method of Lewis and WRISLEY is based on the application 
of a considerable number of shakers. However, according to (2) it is possible 
to excite sufficiently pure natural modes with relatively few shakers. 


(Eingegangen: 26, Februar 1954.) 
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On Stresses of an Isotropic Elastic Disc in the Form of a Cardioid 
Rotating Steadily in its Plane 


By DEBENDRANATH MiITRA, Kharagpur’) 


1. Summary: Two dimensional problems of elasticity have been solved by 
STEVENSON), by expressing the stresses and displacements in terms of two 
complex potential functions. Cases of disc in the form of a circle and of an ellipse © 
rotating steadily have also been solved by him. In this paper the method of 
MUSKHELISHVILI) has been adopted and the equations giving the stress-functions 
have been transformed into integral equations by means of the mapping function 
which represents conformally the region included by the boundary of the disc on 
a circle. The stress functions are then determined by solving the above integral 
equations using function-theory. The case for a cardioid is solved as an example. 


2. A disc is rotating steadily about the origin and the problem is treated as 
one of plain strain. Let D, denote the displacement and 0, ® the stress combina- 
tions be expressed as 


D=u+4+10, (2.1) 
O=4%4+9), (252) 
OD=HxX~VVLQALEY. (2.3) 
From the stress-strain relations it is deduced, 
0D 0D 
(1-29) @=2u (5 + Se), (2.4) 
oD 
O= 4 7 ie 
lu oF (253) 
where 
y) 


the bar placed over a function, denoting the conjugate of the function. 0, ® are 
easily expressed in terms of ’(z), w’(z) which are analytic functions of z only, 


where 0W(z, z)/dz is the body force and y, the elastic constant is given by 
y= 2(1— 2 %)/(1— 7) and dashes over a function indicate derivatives with 
respect to the argument of the function. 


1) Indian Institute of Technology. 
*) A.C. STEVENSON, Phil. Mag. 34, 766-793 (1943). 
3) N. I. MusKuELtsuvitt, ZAMM. 18, 264 (1933). 
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Denoting the directions of the tangent and the normal at a point on the 
bounding curve of the disc by , s respectively we get using (2.6) and (2.7) 


; * 2 
0 ; a 
0 Te) Sm P(e) 4 DO + we ye). (2.8) 


4(tm+im 


If the angular velocity of the rotating disc be 1, 


; 2» 5 Ot pee. 
OW es ae - Ws 


Oz 2 (2.9) 


Also since the boundary of the rotating disc is unstressed, 7” = 53 = 0, and from 
(2.8) the boundary condition in the form given by STEVENsoN?) is obtained 


Qe) + 2B + (2) -—yveW=2er [zz dz. (2. 10) 


3. We now adopt MUSKHELISHVILI’s?) method. Let the mapping function which 
represents conformally the region bounded by the boundary of the disc on the 


- unit circle be 


NA ee ee Te ee ee ee 
’ hase a Le] © v 


z= m(C). (3.1) 


Putting Q(z) = Q{m(C)} = Q(2) and ¢ =1/f on the unit circle, from (2.10) we 
have on the unit circle, 


arty + eke) oA =yoW(t,+20r . m(t) m2) m’(e) ac! 


m’(é) m'(E (3. 2) 


eit (GC) MESON. 


The functions 2(¢), w() can be determined by evaluating the two following 
integral equations obtained by multiplying (3.2) and its conjugate by 
1 dC 
2%% C—«o 


and integrating on the circumference, c of the unit circle. 


gat [20 (e) aay + rae] Ce) ta 


c ce 


zai | *(z) (¢ = = zat] FO m (=) an 


zerf OF) Oat sar fe) LO Goa 


c 


Pa aera tat! F(z) m0) Goa 


Cc 


a 
—~ 
Ww 
aS 
Za 


where o is any point within the unit circle. 


1) A.C, SreveNSON, Phil. Mag. 34, 766-793 (1943). 
2) N. I. Muskuecisuvit1, ZAMM. 13, 264 (1988). 
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4. As an example, we work out the case for a steadily rotating disc in the © 
form of a cardioid. The formula transforming the cardioid on the unit circle is — 


z=m(C)= —a(1 — fa \ es (45) ; 
FQ) =—Lynart zon [225 d= pyena Sa 
_(4. 2) 


a 
—4ora3 (Fo - ZF et4 58-104 5 + Slogé), 


c 
perellosiiigies Wrats 


c 


= —29rtat[—1 (ot— 70% + 210%— 356 + 35) (4. 3) 
1 23 20 (1) | 
4 3 2 a 23 Pa st ! 
+4 {io ie oan ZS 150+ 10—5 = log(t + ot], 
1 yell re : 
es [@ ( Bas ac, [the constant term in w’(¢)]. (4.4) 
On the assumption of the {2(¢) in the form of a power series, 
OG) Sa 4a, 62 + a, 0? (4.5) 
when the constants are to be determined, 
1 ? dt af =,(1 an, ae 
saz [0 m0 + a [ O(E) me eS 
c . és (4.5) 
= —24a Q(o) ——- a{a, (o — 1)? + 2a, (0 — 2) + 3 a3}. | 
The equation (3.3) with (4.3)—(4.6) gives 
o { 2a,—34a 65 o 
Q(o) = —-— 1 2 2 = 2 
i Ft ene ike 21 0) 
{— (1/4) y (o& — 704+ 21 68 — 3567+ 350) 
+ or? a3 ———+-* ee Re = — aa eset) 
pyc CP!) Cu ae bia ie hee CR RS ERE ely 
(1 — o) r. aah ao ie 
Since 2’(0) = a,, Q"(0) = ay, 2(0) = a3, we have from (4.7), 
3@,;—44,+34,=¢,+ 207% a? (- = y + 20), | 
a, + 4g— 3a, = —C, + 20 0 72 a? , (4. 8) 
6 ay — 10 ag= —3 0+ 9 grtaty, | 


The equation (3.4) is evaluated and the following results obtained, 


1 (1 a 
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Be fot n(2) A, ~ [taanee netae 28) aus 
c 

; raz f oO eos Sea); (4.11) 
c 

F s (4.12) 


— 4972 a? {o? —o — 545 (1— 9) log(1 + 9)}| 


(3.4) with (4.9) — (4.12) gives 


; 1 = 0)2,Q’ 2 
a'(o) = 2aa,+al' - Cm ayes 4 — 


—2agreat|% (o®— 70% + 210 ~ 35) (a5 83) 
— 40% +40 + 20— 20 (1 — 0) log(1 + 0)]. 


_ Putting o = 0 in (4.13) we get 


i 35 
3 dy — 40, + 30, = Cot 2er2a*(—— y+ 20). (4. 14) 


The values of the constants are now obtained from (4.8) and (4. 14) 


— 240 8 
f= 0, a, = gra ( ieee 
200 + 203 120 + 231 eeu 
a, = grta® (— 5 ry, a,= ora? (= as). 


Hence (4.7), (4.13) together with (4.15) give Q2(c), w’(o). The hoop stress on the 
boundary given by 


is then easily calculated from (4.7). 
The author has also worked out the case of an inverse of an ellipse by the: 


- above method. 
Zusammenfassung 


Die zweidimensionalen Elastizitatsprobleme wurden von STEVENSON (1943) 
geldst, indem er die Spannungen und Verschiebungen durch zwei komplexe Po- 
tentialfunktionen ausdriickte. In der vorliegenden Arbeit wurde die Methode von 
MuSKHELISHVILI angewendet, um die Gleichungen, welche die Spannungsfunk- 
tionen liefern, als gleichwertige Integralgleichungen darzustellen. Hierzu muss 
die Formel, welche den vom Umfang einer Scheibe eingeschlossenen Bereich auf 
einem Kreis konform abbildet, bekannt sein. Die Spannungsfunktionen sind dann 
durch die Lésung der obigen Integralgleichungen mittels der Funktionentheorie 
festgelegt. Der Fall des Kardioids ist als Beispiel gelost. 


(Received: September 29, 1954.) 
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Hydrodynamisches Modell zur Darstellung von Diffusionsvorgangen 
mit konzentrationsabhangigen Diffusionskoeffizienten 


] 
: 
4 
é 


Von WERNER KNAPPE, Darmstadt?) 


Bei der mathematischen Behandlung von Diffusionsvorgangen geht man im 
allgemeinen vom I. Fickschen Gesetz aus: 


S = —D gradc, (1) 


S = Diffusionsstromdichte (gcm7? s7}), 
c = Konzentration (gcm7), 
D = Diffusionskoeffizient (cm? s—}), 


Hieraus folgt durch Bildung von div G. 


= —divS = div G(D grado), (2) 


ip IKE (US) 


Nur im Sonderfall eines konstanten Diffusionskoeffizienten ergibt sich das soge- 
nannte II. Ficksche Gesetz 


oo 
in 


Wahrend zu (3) zahlreiche Lésungen bekannt sind, gestaltet sich die Anwendung 
der nichtlinearen partiellen D-Gleichung (2) schwierig. Fiir den Sonderfall eines 
beiderseitig unendlich ausgedehnten Systems ist eine exakte Lésung nach Bottz- 
MANN*) méglich. Bei anderen Rand- und Anfangsbedingungen ist man auf Nahe- 
rungsverfahren*) angewiesen. Die damit verbundene erhebliche Rechenarbeit 
legt es nahe, nach physikalischen Modellen zu suchen, in welchen dem Diffu- 
sionsvorgang analoge Prozesse ablaufen. Hierbei liefert die der Konzentration 
entsprechende Messgrésse die gesuchten Zahlenwerte. Modelle dieser Art sind 
schon mit Erfolg zur Behandlung von mathematischen Problemen in der Technik 
eingesetzt worden‘). So lassen sich zum Beispiel instationadre Warmeleitungsvor- 
gange mit Hilfe einer Kette von elektrischen Widerstanden und Kondensatoren 
darstellen, wobei die Spannung der Temperatur entspricht. Ein im Aufbau ahn- 
liches Modell auf hydrodynamischer Grundlage ist von A. D. Moore’) angegeben 
und von A, Doprowsky zum Lésen axialsymmetrischer instationarer Diffusions- 
probleme benutzt worden’). 

Zylindrische Behalter Z aus Glas (Querschnitt = F) werden durch Kapillaren K 
zu einer Kette miteinander verbunden (siehe Figur 1). In einem solchen System 
entspricht der Fliissigkeitsstand h der Konzentration. Die Analogie zum I. Fick- 
schen Gesetz ist gegeben durch das Hagen-Poiseullesche Gesetz, welches den Zu- 


D div gradc = D Ac. (3) 


1) Deutsches Kunststoff-Institut an der Technischen Hochschule. 
L. Botrzmann, Ann. Phys. Chem. 53, 959 (1894) 


i 


) 

) J. CRANK und M. E. Henry, Trans. Far. Soc. 45, 636 (1949). 
) F. Bruckmayer, Allg. Warmetechn. 4, 79 (1953) u. a. 
) 
) 


oo fF Oo 


A. D. Moors, Ind. Eng. Chem. 28, 704 (1936) 
A. Dosrowsky, Kolloidz. 126, 1 (1952). 


J 


a 


. 


_ 
‘ 


. 
¥ 
' 


4 
5 
, 
: 
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bzw. Abfluss zu den einzelnen Behaltern regelt 


5 - Apart _ (haar Bn) OB m1" 
alia > ee rice SC: (4) 


Q,, = Flissigkeitsvolumen, welches pro Zeiteinheit in den n-ten Behal- 
ter Z einstr6mt (cm’s7}), 
p = hydrostatischer Druck (gs~* cm7}), 
vy = Radius der Kapillaren K(cm), 
n = Viskositat der strémenden Fliissigkeit (gs~t cm™?). 


Hila ae ee 
ey ates 
UAT at a aac 


' 


Pic din 


Figur 1 
Skizze des hydrodynamischen Modells fiir konstanten Diffusionskoeffizienten. 


Die Anderung des Fliissigkeitsstandes / in den einzelnen Behdltern ergibt sich zu 


Se See rig ae LS ~ On 
Ot oa F (Qn, Qn-1) cin ps { Ah, — Ahn1} a wae ; Ax? ag Ox2 ° (5) 


Gleichung (5) stellt also das II. Ficksche Gesetz in Differenzenform fiir den 


- inearen Fall dar. Das Modell liefert demnach Naherungslésungen zu ey: 


Ersetzt man die zylindrischen Behalter durch solche, deren Querschnitt F 


von der Hohe h abhangt F = F(h) (siehe Figur 2), so geht (5) tiber in 


iia 2h Aes ohn 
ot «=F(h) | «Ax? Ox? * 


Das Modell arbeitet also mit einem von der Hohe h des Fliissigkeitsstandes und 


D*(h) (6) 


_ damit von der Konzentration abhingigen Diffusionskoeffizienten. 


Figur 2 
Skizze des hydrodynamischen Modells fiir konzentrationsabhangigen Diffusionskoeffizienten. 
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Fiir den linearen Fall erhalt man aus Gleichung (2) allerdings einen von (6) 
abweichenden Ausdruck, namlich 
0c 
ap TA pps, eee 
Es lasst sich aber zeigen, dass bei kleinen Konzentrationsunterschieden (Uber- 
gang zur Differenzenform und Ac > 0) das 2. Glied auf der rechten Seite rascher 
gegen 0 geht als das erste. Wir k6nnen schreiben 


0c 0D (se) 
Nagy 3" 


0c 02c a ( 0c 
fie |= 
Ox 


)>o. (8) 


Figur 3 


Aufbauskizze fuir eine praktische Ausfiihrung. 


Die Annaherung an (2) wird also mit fortschreitender Unterteilung der Kette und 
damit kleineren Hohendifferenzen zwischen den Fliissigkeitsspiegeln der einzelnen 
Behalter besser, Unerwiinscht hohe Differenzen treten zwangslaufig zu Beginn 
eines Diffusionsvorgangs auf. Um Fehler zu vermeiden, verfahrt man am besten 
so, dass man die Konzentrationsverteilung fiir das erste Stadium des Diffusions- 
vor gangs nach der Methode von BoLttzMAnn}) berechnet und die exakten Werte 
als Anfangswerte fiir das Modell verwendet. 

Ein Modell in der angegebenen Art lasst sich im Prinzip folgendermassen her- 
stellen (Figur 3): Man unterteilt einen rechteckigen Trog durch Einbringen zahl- 
reicher senkrecht stehender Platten P in viele gleich grosse schmale Zellen. 
Durch Einlegen einer entsprechend F = F(h) geschnittenen Folie G in jede dieser 
Zellen lasst sich das Modell dem jeweiligen Problem anpassen. 


Summary 


It is possible, to get approximate solutions of the differential equation for 
non-stationary diffusion (F1cKs 2nd law) by means of a hydrodynamic model, 
consisting of cylindric containers connected by capillary tubes. Solutions with 
a diffusion coefficient dependent of concentration may be obtained by taking 


containers with a cross section variable with the height. Some proposals for the 
design of such a model are made. 


(Eingegangen: 30. September 1954.) 
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Varia — Miscellaneous — Divers 


Eine Bemerkung zur Mitteilung von H. Wundt: 
«Uber eine Funktionalgleichung aus der Warmeleitung’) » 


Von Mixtés Hosszu in Miskolc?) (Ungarn) 


Herr H. Wunpt behauptet in der oben stehenden Arbeit, dass die Funktion 


x 


f(x) Aft +B (7) 


t— Int—1 
1 


eine Lésung der Funktionalgleichung 


(ee ee ge (1) 


(7) die Gleichung (1) auch wirklich 


ist. Er hat aber nicht naher untersucht, ob 
befriedigt. 

Wir werden zeigen, dass die Funktio 
Lésung als die triviale f(¥) =¢ (const) h 
der Gleichung (1) nach ¥ bzw. y: 


(e) = Af(2) baw. y= PO): 


nalgleichung (1) keine andere derivierbare 
at. Wir bilden die partiellen Ableitungen 


Dies gibt eine notwendige Bedingung*) fiir die Existenz einer nichttrivialen 


Lésung der Funktionalgleichung (iis 

0? ae e 
Ox Oy = Ly Ox Oy 
)/(In #/y) ist, wird diese Bedingung nicht erfiillt, Um dieses 


D(x, y) = pou ey tea 7"(9)) 1 10" 


Wenn z= (*—-— Vy 


zu zeigen, bilden wir 
z, _ Ina/y—1/x (¥— 9) 
zy —ima/y+ify(#—y)- 
und weiter 
0 in Z# ere yee Tat Lives 
(BR 8 —Inz/y—1+4+ y/% —Ins/y+4/y—1’ 


und 
—1/x? (Inx/y —4-+ y/*) — (A/a — y/¥") (—1f/y + 1/¥) 


ee) (inz/y -1+9/*)? 
—1/y? (—In#/y + #/y — 1) — (—1/% + 1/y) (A/y — #1") 
t (—Ina/y + /y — 1)? 


1) Erschienen in der ZAMP 4, 172-175 (1954). 
2) Mathematisches Institut der Technischen Universitat. 
3) P, Saint ROBERT, De la résolution de certaines équations & trois variables par le moyen dune 


végle glissante, Mem. R. Accad. Torino (2) 25, 538-72 (1871). 
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Setzt man zum Beispiel 7 = 1, y = 1/e, dann wird 


_ —(1-1 —etl (= ed) = [BS ee ee 
ION Ny YB) zd A ane — ; ae ie res 
(e — 2)*[—e — e (e —1) (1—e)] + ele? — 2e + (e —1) (1 — @)] 
— jae (e — 2)? 
graye lle= Ae 24) =e ele — 2 i 
e 
Teaaye P,(e) +0, 


da e transzendent ist, also keine algebraische Gleichung [wie P,(e) = 0] befrie- 
digen kann. 


Da die Konstante f(x) = ¢ (1) befriedigt, ist sie die einzige differenzierbare _ 


Losung von (1). 


Summary 


This paper shows, in connection with H. Wunpt’s paper: Uber eine Funk-— 


tionalglerchung aus der Warmeleitung [ZAMP 5, 172-175 (1954)], that the only 
differentiable solution of the functional equation 


fie Bg ye 1 
2 (Rae) = He) + Ho) (1) 
is the trivial /(¥) = constant, and thus H. Wunpv’s solution 


x 


Int+ 1/t—1 
I(x) a [SS ih 


1 


with A + 0 does not satisfy (1). 


Bemerkung zur obigen Mitteilung von Herrn M. Hosszu 


Herr Hosszu hat mit Recht festgestellt, dass die von H. Wunpt vorgelegte 
Funktion 
Se NER aici 
In’ — Iny (1) 


nicht durch ein Leiternomogramm mit drei geradlinigen Leitern dargestellt wer- 
den kann. Herr Hosszvu zeigte zu diesem Zweck, dass (1) die Saint-Robertsche 
Differentialgleichung nicht erfiillt, was man natiirlich auch einsehen kann, ohne 
die Transzendenz von e zu bemiihen. Es mag von Interesse sein, anzugeben, wie 
das Leiternomogramm von (1) wirklich aussieht. Bezeichnet man mit &, 9 recht- 
winklige Koordinaten in der Ebene, so lauten die drei Skalen des Nomogramms 


nf x 


x-Skala 6 = or , = ina : (2) 
I y 

/-Skal 5 cay Se aye 

4 ala € ing? n gh (3) 


z-Skala €&=0, n= 2). (4) 


"a 


ite) See —_— ROT TRA Sy eee a ee ee 
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‘Man bemerkt, dass die ersten beiden Skalen zusammenfallen und dass die z-Skala 
geradlinig ist. Die Bedingung, dass die drei Punkte (2) (3) (4) auf einer Geraden 
liegen, ist in der Tat identisch mit der vorgegebenen Beziehung (1). 

Die Behauptung von Herrn Hosszuv lasst sich nun ohne jede Rechnung aus 
folgendem Satz der Gewebetheorie herleiten: 

«Ein Leiternomogramm kann nur dann in ein Nomogramm mit drei geraden 
Leitern verwandelt werden, wenn die drei Leitern Stiicke ein- und derselben 
(eventuell zerfallenden) Kurve 3. Grades sind.» Da die Skala (2) transzendent ist, 


kann also (1) nicht durch drei geradlinige Leitern nomographiert werden. 
E. STIEFEL 


(Eingegangen: 19. Oktober 1954.) 


Kleine Bemerkung zur asymptotischen Entwicklung 
des Fehlerintegrals 


Von PETER HENRIcI, Washington, D.C.*) 


In einer kiirzlich in dieser Zeitschrift erschienenen Mitteilung [6]?) gibt 
J. ZBoRNIK eine asymptotische Entwicklung fiir ein Integral, das auf die Form 


fo, 2) 
ole) = fete ar (jarge] S29 
a 
gebracht werden kann und daher nahe mit dem Komplement des Fehlerintegrals 


verwandt ist®). Im Hinblick auf einige Ungenauigkeiten in [6] sind vielleicht 


einige erganzende Angaben am Platze. 
1. Durch einfache Variabelntransformation, Entwicklung des Binoms und 
gliedweise Integration findet man unmittelbar das bekannte, mit demjenigen von 


[6] identische Resultat 


ey = few (at g) Maem ae 
(5), i a (> + 1) ur (J n= 1) 


gesetzt ist. Der asymptotische Charakter der Reihe folgt dabei zum Beispiel aus 
dem bekannten Watsonschen Lemma (vgl. [5], S. 236), nicht aber natiirlich wie 


- pei [6] durch den blossen Nachweis, dass das Restglied fiir z > co gegen Null 


strebt. 
2. Auch der in [6] der Fehlerabschatzung zugrunde gelegte Satz («In asympto- 


tischen Reihen kann der Betrag des Restgliedes den des letzten Gliedes der Reihe 


1) American University. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 146. 


3) Das Integral war bei einer Untersuchung von H. Hésri tiber den Einfluss des Windes auf 
Raketenbahnen aufgetreten. 


ZAMP VI/10 


146 Varia — Miscellaneous — Divers ZAMP 


nicht iibersteigen») ist unrichtig!). Die numerische Abschatzung des Restes erfor- 
dert in den meisten Fallen eine besondere Untersuchung. Im vorliegenden Falle 
lasst sich, wie hier ohne Beweis kurz mitgeteilt werden soll, folgendes zeigen: Es ist 


e7* | "oy (1/2), (1/2)n | 
Ve ey art tO ear 


p(2) = 


mit 
a —n—(1/2) 
o—fe(1+ 5) a ds; 
0 


man findet, z = v e*® gesetzt, fiir || S (2/2) 
|6| <1, sign(arg@) = sign#, |arg9| S >: 
Die angegebenen Schranken sind, jede fiir sich (aber nicht notwendig beide gleich- 
zeitig), exakt. Um einen méglichst guten Wert fiir y(z) zu erhalten, wird man also 
das letzte beriicksichtigte Glied der asymptotischen Reihe mit (1+ 7 sign #)/2 
multiplizieren; der Fehler ist dann absolut kleiner als das so modifizierte Glied?). 
3. Andere numerisch wertvolle Verfahren zur Berechnung des Fehlerintegrals 
wurden von SCHLOMILCH [4] (Kettenbruch) und LaIBLe [2] angegeben. Eine Tafel 
fiir komplexe Werte von z ist von CLEMMow und MUNDFoRD [1] berechnet worden. 


LITERATURANGABEN 


[1] P. C. CLemmow and C. M. MunpForp, Phil. Trans. Roy. Soc. London [A] 245, 
189-212 (1952-1953). 

[2] T. Larsie, Z. angew. Math. Physik 2, 484—486 (1951). 

[3] F. Léscu und F. ScHosiik, Die Fakultadt (Leipzig 1951). 

[4] O. ScuLt6mitcnH, Z. Math. Phys. 76, 261—262 (1871). 

[5] G. N. Watson, Bessel Functions, 2nd ed. (Cambridge 1944). 

[6] J. ZBoRNIK, Z. angew. Math. Phys. 5, 345—351 (1954). 


(Eingegangen: 6. Januar 1955.) 


Internationales Symposium iiber elektrische Entladungen in Gasen 


Vom 25, bis 30. April 1955 wird in Delft (Technische Hochschule), Nieder- 
lande, ein internationales Symposium iiber elektrische Entladungen in Gasen 
stattfinden. Die Mittel fiir die Durchfiihrung dieser Veranstaltung werden 
einerseits im Rahmen der Unesco von der Union Internationale de Physique 
Pure et Appliquée, anderseits von der Technischen Hochschule Delft und von 
den Philips-Werken, Eindhoven, aufgebracht. 

Die folgenden Hauptthemen sollen diskutiert werden: 


1, Fundamental-Prozesse; neue Ansichten tiber die Gasentladungsmecha- 
nismen. 


1) Ein besonders drastisches Gegenbeispiel ist die identisch verschwindende asymptotische 
Entwicklung von e~* fiir x -> oo, 

2) Das hieraus fiir reelle Werte von z (# = 0) ergebende Resultat ist wohlbekannt, Fiir kom- 
plexe z leistet unsere Abschatzung im betrachteten Winkelbereich mehr als diejenige in [3], S. 129f. 


ae ee 
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. Instabilitaten, Stabilitatsbedingungen und Plasmaschwingungen. 
. Ziindspannung als Funktion von p.d. und Frequenz. 
. a) Moderne Messmethoden, angewendet auf Gasentladungsuntersuchungen. 
b) Anwendung der Gasentladung zu Messzwecken in anderen Gebieten. 
. Bogenentladungen. 
Funken. 
. Verschiedenes (Geiger-Miiller-Réhren, Ionenquellen). 
ereits habensich die folgenden Referenten zu einem Hauptvortrag bereit erklart: 
Prof. H.S.W.Massey, London: «Fundamental Primary Processes in Gas 
Discharges. » 
Prof. L.B.Lors, Berkeley: « Field Measurements in Glow Discharges with 
‘a Refined Electron Beam Probe and Automatic Recording.» 
Ir CHR. VAN GEEL, Delft: « On the Internal Selfinduction of Gas Discharges.» 
Prof. S.C. Brown, Cambridge, Mass.: « Breakdown in Gases at Microwave 
Frequencies. » 
Dr. M.A. Bronp1, East Pittsburgh Pa.: « Microwave and Optical Techniques 
for Gas Discharges.» 
Prof. W. LocutE-HoLtGREVEN, Kiel: «Bogenentladungen.» 
Prof. J.M.MeeEk, Liverpool: « Spark Discharges.» 
Die zugelassenen Sprachen des Symposiums sind Franzésisch, Deutsch und 
Englisch; in der Diskussion ist Englisch bevorzugt. 
Zwei Tage sind reserviert worden fiir Exkursionen. Fiir die Damen der 
Symposiumteilnehmer ist ein spezielles Programm vorgesehen. 
Fiir Auskiinfte iiber das Symposium wende man sich an den Sekretar des 
Symposiumausschusses Ir A. W.VAN WAGENSVELD, Mijnbouwplein 11, Delft, 


Niederlande. 
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Gas Turbine Analysis and Practice. By BurGEss H. JENNINGS and 
WitLaRD L. Rocers (McGraw-Hill Publ. Co., Ltd., London, 1953). 475 pp., 
200 Fig.; 64s. 

Das Buch gibt zunachst einen kurzen geschichtlichen Riickblick auf die Ent- 
wicklung der Gasturbine und geht gleichzeitig auf ihre Wirkungsweise ein. Es 
folgen Abschnitte tiber die thermodynamischen Grundlagen im allgemeinen (ins- 
besondere iiber die Eigenschaften der Luft), die Prozessfiihrung der Gasturbine 
und die Theorie der Turbinen und Verdichter. Im beschreibenden Teil sind sowohl 
Flugtriebwerke als auch Gasturbinen fiir andere Zwecke behandelt. Es folgen 
weitere Kapitel iiber Verbrennung und Brennkammer, Festigkeitsprobleme, 
Schwingungsprobleme und metallurgische Fragen. 

Der allgemeine Charakter des Buches ist elementar. Trotzdem kénnte die 
thermodynamische Einleitung mehr in die Tiefe gehen. Unmodern und geradezu 
primitiv ist die Behandlung der Turbinen- und Verdichtertheorie. Der beschrei- 
bende Teil ist sowohl fiir die Flugtriebwerke als auch fiir die stationaren und 
Fahrzeuggasturbinen knapp gefasst,. wogegen der Problemkomplex «Brennstoff, 
Verbrennung, Brennkammer» eine fiir ein elementares Buch recht ausfiihrliche 
Behandlung erfahrt. Allerdings' lasst die Art:der Betrachtung dieser Fragen er- 
kennen, dass, vom Standpunkt des Flugtriebwerkbauers ausgegangen wird. So 
kommt das heute fiir den tibrigen Gasturbinenbau so wichtige Problem der 
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Schwerélverbrennung, das mit den Stichworten Vanadium, Natrium, Schwefel 
angedeutet werden kann, vollstandig zu kurz. 


Die heute bekannten Lehrbiicher des Gasturbinenbaues befriedigen zumeist | 


nicht, und dies gilt auch von demjenigen von JENNINGS und Rocers. Man muss 


allerdings beriicksichtigen, dass der amerikanische Autor von anderen Gegeben- 


heiten ausgehen muss als etwa der kontinentaleuropaische. Er muss Spezialkennt- 
nisse vermitteln, ohne die gleiche allgemeine Schulung voraussetzen zu k6nnen 
wie bei uns. Dies ist dem vorliegenden Buch deutlich anzumerken. Auf eine nicht 
zu unterschaétzende Qualitat sei iibrigens noch hingewiesen. Das Buch enthalt 
namentlich in den Abschnitten tiber Thermodynamik und Verbrennung eine 
Fiille wertvoller Daten in Tabellen zusammengefasst (im englischen MaSsystem), 
die man sonst unter Umstanden miihsam zusammensuchen misste. W. Traupel 


Lexikon der Physik. Herausgegeben von H. FRANKE (Franckh’sche Ver- 
lagshandlung, Stuttgart, 1950). 2 Bande, 720 und 828 S., mit 1263 graphischen 
Darstellungen und 232 Tafelbildern; I. Bd. DM 72.-, II. Bd. DM 82.-. 

Es ist leider richtig, wenn im Vorwort des in angemessener Aufmachung er- 
schienenen, in zwei Bande aufgeteilten Lexikons der Physik vermerkt wird, dass 
die betrachtliche Ausweitung, die das Gebiet der Physik erfahren hat, notgedrun- 
gen zu einer starken Spezialisierung fiihren musste. Damit hat sich fiir Physiker 
und Ingenieure das Bediirfnis nach einem Nachschlagewerk eingestellt, das ihnen 
erlaubt, sich schnell und trotzdem zuverlassig tiber Belange benachbarter Wis- 
sensgebiete zu informieren. Diesem Ziele kommt das vorliegende Lexikon in er- 
freulicher Weise nahe. 

Das Werk behandelt rund 14000 Gegenstande. Auswahl, Bezeichnung und 
Umfang der Erérterungen sind unvermeidlicherweise eine Ermessensfrage. Jeden- 
falls war es aber richtig, grundsatzlich auf eine Schilderung der Messverfahren 
und Versuchseinrichtungen zu verzichten, um so den Umfang des Werkes in er- 
traglichen Grenzen zu halten; die Hinweise auf die einschlagige Literatur geniigen 
und sind ausserst wertvoll. 

Unter den Gegenstanden finden sich auch Angaben iiber das Leben der mar- 
kantesten Forscherpersénlichkeiten der Physik. Der Rezensent vermisst darunter 
Namen wie BLocu, voN KARMAN, insbesondere APPLETON. Doch ist es zugege- 
benermassen schwer, die Grenzen einer verniinftigen Auswahl zu finden. Auch 
sei vermerkt, dass die Begriindung fiir die Erteilung des Nobelpreises an EINSTEIN 
unrichtig ist; ebenso die Aussage, dass zum Beispiel DEBYE zwei Nobelpreise 
(Physik und Chemie) zuerkannt wurden. Diese Beanstandungen sind aber sicher- 
lich von ganz untergeordneter Bedeutung gegeniiber dem grossen Vorteil und 
Dienst, den das Lexikon der Physik jedem Interessenten zu erweisen vermag. 

R. Sanger 


Formeln und Tabellen der mathematischen Statistik. Von U. Grar 
und H. J. Hennine (Verlag Springer, Berlin 1953). 102 S.; DM. 9.—. 

Gut iiberlegte, rezeptmassige Zusammenstellung der wesentlichen Formeln, 
Tabellen und Nomogramme fiir den Gebrauch der wichtigsten statistischen 
Methoden und Teste in der Praxis. Die Formeln sind ohne Beweise dargestellt, 
hingegen werden Beispiele formuliert und diskutiert. Das Biichlein enthalt auch 
die Grundlagen fiir Stichprobenplane und Kontrollkarten. 

Die Formelsammlung ist fiir Forscher und Ingenieure mit den notigen theo- 
retischen Kenntnissen in der mathematischen Statistik ein wertvolles technisches 
Hilfsmittel, numerische Ergebnisse statistisch zu werten und zu analysieren. 


W. Saxer 
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Einfiihrung in die Mikrowellen-Elektronik. Band 16 der « Monographien 
der elektronischen Nachrichtentechnik», Von W. KLEEN. (S. Hirzel Verlag, Ziirich 
#952)300 'S.;,129-Abb.; sFr. 30.—. 

Angesichts der sehr raschen Entwicklung der Hochfrequenztechnik und der 
Elektronik, namentlich im Bereiche héchster Frequenzen, wo sich diese beiden 
Gebiete aufs innigste durchdringen, besteht ein grosses Bediirfnis nach zusam- 
menfassenden Werken, welche kurz, klar und doch mit einer gewissen Vollstan- 
digkeit und Griindlichkeit die jiingsten Fortschritte zur Darstellung bringen. 
Diesem Ziele gerecht zu werden, ist W. KLEEN in seiner Einfihrung in die Mikro- 
wellen-Elektronik aufs schénste gelungen. 

Aus dem sehr weitverzweigten Stoff die richtige Auswahl zu treffen und den- 
selben in einem systematischen Aufbau iibersichtlich zu gliedern, bedeutet eine 
miihsame Aufgabe, die in vorbildlicher Weise gelést ist. Eine kurze Aufzahlung 
der Kapitel mége eine Vorstellung von dem reichen Inhalt vermitteln: 1. Die 
Systeme der Mikrowellen-Elektronik; 2. Elektronenbewegung in statischen Fel- 
dern; 3. Leistungsaustausch zwischen Elektronenstromungen und elektrischen 
Wechselfeldern; 4. Str6me in Mikrowellen-Rohren; 5. Geschwindigkeitsmodulation 
in stehenden Feldern; 6. Phasenfokussierung in von Hochfrequenz freien Rau- 
men als ballistisches Problem; 7. Leistungsauskopplung aus stehenden Feldern; 
8. Raumladungssteuerung; 9. Phasenaussortierung; 10. Steuerung von Elektro- 
nenstrémungen durch fortschreitende Wellen bei Abwesenheit statischer Quer- 
felder; 11. Phasenfokussierung und freie Raumladungswellen; 12. Wechsel- 
wirkung zwischen Elektronenstromungen und fortschreitenden Wellen in ge- 
kreuzten statischen Feldern; 13. Systematik der Mikrowellenréhren; 14. Tech- 
nische Bedeutung und Anwendung von Mikrowellenrohren,; 15. Die R6hre in 
der Schaltung; 16. Rauschen; 17. Schwingungskreise ; 18. Verzégerungsleitungen ; 
19. Elektronenkanonen. 

Der Raum verbietet leider, auf Einzelheiten einzugehen. Die Darstellung ist 
knapp, aber fasslich gehalten. Gewisse mathematische und theoretische Grund- 
lagen, iiber welche der Nachrichteningenieur verfiigen soll, mtissen selbstver- 
standlich vorausgesetzt werden. Kiinftigen Auflagen wiirden wir die Verwendung 
des Giorgi-MaBsystems empfehlen. Hervorgehoben seien die vielen aufschluss- 
reichen Figuren; willkommen ist die jedem Kapitel am Schluss beigegebene Lite- 
raturiibersicht. Die Ausstattung des Buches ist ausgezeichnet. Mége es viele 
Leser finden. F. Tank 


Experimental Nuclear Physics. Vol. I. By Secrt (J. Wiley, Inc., New 
York, 1953). 780 pp.; $15.00. 

Dieses Buch stellt den ersten Band einer Handbuchfolge dar, in welcher eine 
zusammenhangende Darstellung unserer gegenwartigen Kenntnisse auf dem Ge- 
biete der Kernphysik vom experimentellen Standpunkt versucht wird, Der vor- 
liegende erste Band enthalt fiinf Beitrage aus verschiedenen, in sich abgeschlos- 
senen Arbeitsgebieten. Titel und Autoren dieser Monographien sind folgende: 
Nachweismethoden (H. H. Stavs, 165 S.); 

Durchgang von Strahlung durch Materie (H. A. BETHE und J. ASHKIN, LOUIS) 5 
. Kernmomente und Statistik (N. F. Ramsey, 109 S.); 

. Zweikérperprobleme und Elemente der Kernstruktur (N. F. Ramsay, 90 S.); 
Dynamik und Optik geladener Teilchen, Isotopenhaufigkeit und Atommassen 
(K. T. BAINBRIDGE, 207 =) eet: 

1. Der Aufsatz iiber Nachweismethoden behandelt die Grundprinzipien der 
verschiedenen Zahlertypen, soweit sie Gaszahler und Ionisationskammern be- 
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treffen, ferner die zur Verstarkung der primaéren Impulse benétigten elektroni- 
schen Schaltungen. Er gibt einen vorziiglichen, knappen Uberblick iiber den 
gegenwartigen Stand der experimentellen Technik auf diesem Gebiet. Die aller- 
neuste Entwicklung, die vor allem durch das Aufkommen der Szintillationszahler 
und Detektoren fiir Cerenkov-Strahlung gekennzeichnet ist, findet allerdings 
keine Beriicksichtigung mehr, sie soll wahrscheinlich einem der folgenden Bande 
vorbehalten bleiben. 

2. Eine ausgezeichnete Zusammenfassung bietet der Artikel iiber die Wech- 
selwirkung zwischen Strahlung und Materie. Er behandelt vor allem das Problem 
der Energieverluste geladener Teilchen und elektromagnetischer Strahlung, und 
zwar in einer fiir die Praxis des Experimentators sehr geeigneten Form. Zahl- 
reiche Tabellen und Kurven kommen den praktischen Bediirfnissen des auf 
diesem Gebiete Arbeitenden entgegen. 

3. Die Monographie tiber Kernmomente und Statistik fiillt eine Liicke in der 


eaten 


bestehenden Literatur. Sie enthalt unter anderem sehr gut brauchbare Zusam- ~ 


menstellungen aller wichtigen experimentellen Daten auf diesem in rascher Ent- 
wicklung begriffenen Arbeitsgebiet. 

4. Auch der zweite Aufsatz des gleichen Autors, welcher das kernphysikalische 
Zweikérperproblem zum Gegenstand hat, bietet eine sehr klare Ubersicht des 
experimentellen Materials und seiner gegenwartigen theoretischen Erklarungsver- 
suche. Es ist allerdings zu bemerken, dass dieses Thema auch schon in anderen 
zusammenfassenden Darstellungen eine erschépfende Bearbeitung gefunden hat. 

5. Der letzte Artikel (BAINBRIDGE) schildert in souveraner Weise die gesamte 
technische Entwicklung wie auch die Theorie der Massenspektroskopie, von ihren 
Anfangen bis in die jiingste Zeit. Auch diese Abhandlung enthalt zahlreiche, fiir 
den Praktiker unentbehrliche Zusammenstellungen der in der Originalliteratur 
verstreuten experimentellen massenspektroskopischen Daten. 

Dieses Sammelwerk ist seiner ganzen Anlage nach fiir den vorgeriickten, zu 
selbstandiger Forschungsarbeit befahigten Leser bestimmt, es wird einem solchen 
sehr gute Dienste leisten. A. Waffler 


Physikalische Grundlagen der Messtechnik in der Warmewirtschaft. 
Von S. VALENTINER. (Verlag fiir angewandte Wissenschaften, Wiesbaden 1953). 
48 Fig., 150 S.; 2. verdanderte und erweiterte Auflage, Leinen, DM 12.80, broschiert 
DM 8.80. 

In vielen technischen Betrieben besteht die dringende Aufgabe, Energie in 
Form von Warme méglichst vollstandig zu nutzen, um so den Betrieb wirt- 
schaftlich zu gestalten. Dies erfordert zuerst eine genaue Kenntnis der Energie- 
betrage, die bei den vorkommenden thermodynamischen Prozessen umgesetzt 
werden, und insbesondere der Warmemengen, die dabei ungenutzt verlorengehen. 
Hierfiir ist meist die Durchfiihrung zuverlassiger Messungen notwendig, deren 
physikalische Grundlagen zuvor bekannt sein sollten. 

Die physikalischen Grundlagen warmetechnischer Messungen in gedrangter, 
fiir den Ingenieur leicht verstandlicher:- Form zu‘ vermitteln, ist der Zweck des 
vorliegenden Bandchens. Es ist aus Vorlesungen des Verfassers an der Bergaka- 
demie Clausthal hervorgegangen. Dementsprechend ist die Auswahl des behan- 
delten Stoffes besonders im: Hinblick auf,den Hiittenbetrieb getroffen worden, 
dem auch die angefiihrten Beispiele zur Hauptsache entstammen. Der Inhalt des 
Buches ist in fiinf Abschnitte gegliedert; ihre Uberschriften lauten: Temperatur- 
messung; spezifische Warme und Warmeiibergang; Das Verhalten der Gase und 
Dampfe; Str6men von Gasen und Dampfen, Mengenmessung; Feuerungstechnik. 
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4 Im ersten Abschnitt wird vor allem die Messung hoher Temperaturen behan- 
~ delt. Mit Recht werden die physikalischen Grundlagen der verschiedenen Arten 
‘ von Strahlungspyrometern besonders ausfiihrlich dargelegt. Dabei wird auch auf 
die haufig entstehenden Messfehler und die zu ihrer Verhiitung notwendigen 
- Massnahmen hingewiesen. Im zweiten Kapitel werden die Gesetzmassigkeiten 
_des Warmetiberganges leider allzu kurz behandelt. Die Warmeiibertragung ist 
_ doch fiir alle warmetechnischen Probleme von so zentraler Bedeutung, dass eine 
) griindliche Kenntnis ihrer Grundgesetze unerlasslich scheint. Sehr summarisch 
werden in den beiden folgenden Abschnitten auch die allgemeinen thermodyna- 
_ mischen Eigenschaften von Gasen und Dampfen erortert; eine Ausnahme bildet 
die etwas umstandliche Behandlung feuchter Gase. Das letzte Kapitel tber 
_ Feuerungstechnik dagegen vermag wieder eher zu befriedigen; es bringt inter- 
essante Hinweise tiber den Ablauf der Verbrennung in technischen Ofen. 

Durch das Reproduktionsverfahren bedingt, konnten in der zweiten Auflage 
des Bandchens zwar Anderungen gegeniiber der ersten im Text beriicksichtigt 
werden; alle Neuzusatze dagegen mussten in einen Anhang verwiesen werden, 
~ was die Ubersichtlichkeit der Darstellung beeintrachtigt. Druck und Ausstattung 
_des Buches sind im tibrigen vorzitiglich. 

Wenn auch noch manche Wiinsche beziiglich Auswahl des Stoffes und tiefer- 
gehender Behandlung wichtiger Teilgebiete offenbleiben, so darf doch erwartet 
werden, dass das Buch den in der Warmewirtschaft tatigen Ingenieuren von er- 
heblichem Nutzen sein wird. K, Elser 
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Die konfluente hypergeometrische Funktion. Von H. BUCHHOLZ 
(Springer-Verlag, Berlin 1953). 234 Se op Abba DMis6—-: 

‘Das vorliegende Buch ist erschienen als Band 2 der von F. LéscH herausge- 
gebenen Reihe «Ergebnisse der angewandten Mathematik». Es ist sehr zu begris- 
sen, dass H. BucHHoLz es unternommen hat, die wichtigsten Eigenschaften der 
konfluenten hypergeometrischen Funktion darzustellen. Die Funktion tritt bei 
zahlreichen Problemen der Physik und der Mechanik, insbesondere auch bei der 
durch Separation erzielbaren Losung der Wellengleichung in parabolischen Koor- 
dinaten auf, und, was ebenso wichtig ist, sie enthalt als Spezialfalle eine grosse 
Zahl von haufig vorkommenden transzendenten und elementaren Funktionen. 
Sie ist definiert als Lésung der Whittakerschen Differentialgleichung zweiter 
Ordnung, die in bekannter Weise aus der Gau®schen Differentialgleichung 
durch Konfluenz der Nullstellen z = 1 und z = co entsteht und neben der Varia- 
blen noch zwei Parameter aufweist und mit der Kummerschen Funktion 
1F,(«; 8; 2) zusammenhangt. Von der urspriinglichen Whittakerschen Definition 
weicht BucuHorz ab durch die in vielen Fallen zweckmassige Ersetzung des 
hintern Parameters m durch p/2 und durch Hinzufiigung des Faktors 1/I(1 + 4), 
was einmal zu einer einfachern Schreibweise fiihrt und die Funktion ausserdem 
zu einer ganzen transzendenten Funktion auch im Parameter « macht. 

Die Kenntnis der konfluenten hypergeometrischen (oder parabolischen) Funk- 
tion ist in den vergangenen Jahren besonders durch Arbziten von A. ERDELYI, 
C. S. MEIJER, F. TRIcomi und andern entscheidend geférdert worden, so dass eine 
Zusammenfassung der wichtigsten Ergebnisse durchaus wiinschenswert geworden 
ist. Uber die dargestellten Beziehungen gibt wohl am besten eine Inhaltsiibersicht 
Auskunft, und diese lasst gleichzeitig erkennen, was im Buche nicht behandelt 
wurde, Die einzelnen Abschnitte sind: I. Die Differentialgleichung der konfluen- 
ten hypergeometrischen Funktion in ihren verschiedenen Formen und die Defi- 
nitionen der sie lésenden Funktionen. II. Allgemeine Integraldarstellungen fur 
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die parabolischen Funktionen selbst und ihre Produkte. III. Die Asymptotik deg | 
parabolischen Funktionen. IV. Unbestimmte und bestimmte Integrale mit para- 
bolischen Funktionen und einige unendliche Reihen. V. Die den parabolischen 
Funktionen zugehérenden Polynome und unendliche Reihen mit diesen Poly- 
nomen. VI. Die Parameterintegrale in den Beziehungen fiir die verschiedenen 
Wellentypen der mathematischen Physik in den parabolischen Koordinaten. 
VII. Nullstellen und Eigenwerte. | 
In einem ersten Anhang sind die Sonderfalle der konfluenten hypergeometri- 
schen Funktion zusammengestellt und gleichzeitig die wichtigsten Vertafelungen | 
angegeben. Der zweite Anhang besteht in einem sehr umfassenden, sorgfaltig zu- 
sammengestellten Literaturverzeichnis, das nicht nur die rein mathematischen 
Arbeiten, sondern auch die Untersuchungen iiber physikalische, technische und 
statistische Probleme enthalt, die auf die parabolische Funktion fiihren. 
_ Im ganzen betrachtet, handelt es sich um ein sehr willkommenes Werk, das 
vor allem dem Physiker und Ingenieur, aber auch dem Mathematiker grosse 
Dienste leisten wird. Der Druck ist hervorragend, und — was angesichts des 
Formelreichtums nicht so selbstverstandlich ist — die Zahl der Druckfehler ist 
minimal. E. Roth-Desmeules 


Vorlesungen uber die Theorie der Integralgleichungen. Von I. G. PrE- 
TROWSKI (Physica-Verlag, Wiirzburg 1953). 100 S.; DM 7.80. 

Le lecteur trouvera dans cet ouvrage une exposition concise, mais trés claire, 
de la théorie «classique» des équations intégrales linéaires; il n’y trouvera aucune 
application aux problémes aux limites des équations différentielles. 

Le chapitre I donne d’abord quelques exemples de problémes qui conduisent a 
des équations intégrales ; il montre ensuite les analogies entre équations intégrales 
linéaires et systemes d’équations algébriques linéaires puis démontre les théorémes 


de FREDHOLM, d’abord dans le cas des noyaux dégénérés ¥ f,(2) g,(g), puis dans: 
1 


celui des noyaux uniformément continus par un passage a la limite. Le chapitre 
se termine par l’étude des noyaux du type K(p, q)/pq* qui interviennent dans les 
applications les plus importantes (fonction de GREEN). 

Le chapitre II traite des équations intégrales de VOLTERRA. 

Le chapitre III expose la théorie d’E. Scumipr des équations a noyau symé- 
trique et les développements en séries orthogonales qui en découlent. A y signaler 
une démonstration intéressante et peu connue de I. M. GELFAND de l’existence 
d’une valeur fondamentale d’un noyau symétrique. 

Deux annexes traitent de la réduction d’une forme quadratique a la forme 
canonique par une transformation orthogonale et de la théorie des équations 
intégrales 4 noyau de carré intégrable au sens de LEBESGUE, M. Plancherel 


Raumkinematische Grundlagen. Von Rupotr BEYER und ERNST SCHOR- 
NER (Johann Ambrosius Barth, Miinchen 1953). 104 S., 33 Abb.: DM 9.60. 

Das Buch stellt eine Einfiihrung in die raumliche Kinematik dar und ent- 
wickelt in anschaulicher Weise die Grundlagen der raumlichen Getriebeanalyse. 
Die Behandlung ist rein graphisch und durch den bemerkenswerten Versuch ge- 
kennzeichnet, dem Leser die Raumkonstruktionen durch stereoskopische Bilder 
nahezubringen. Die einzelnen Abschnitte sind: 1. Vektorielle und zeichnerische 
Grundlagen der Raummechanik; 2. Zusammensetzung raumlicher Bewegungen; 
Grundlagen fiir Analyse und Synthese der Raumgetriebe, 3. Beschleunigungs- 
verhaltnisse der Schraubenbewegung. H. Ziegler 


Neuerscheinung 
Der Briefwechsel 
von Johann Bernoulli 
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4 Herausgegeben von der Naturforschenden Gesellschaft in Basel ; 
2 Band 1 

3 532 Seiten, 162 Briefe im Originaltext mit 128 faksimilierten Zeichnungen, 

a 3 Schriftproben und 2 Portrits (1955). In Leinen gebunden Fr. 60.-. 

E Die Ausgabe kann auf Verlangen auch broschiett geliefert wetden. 

: 
‘ c 
q Das Vorwort des Herausgebets (O. Spiess) erzahlt die Vorgeschichte des Unter- 
‘ nehmens, insbesondete die Entdeckung des lange verschollenen ‘Bernoullischen 

; Nachlasses. Daran schliesst sich eine Ubersicht tiber das vothandene Material und 

; eine Skizze des Editionsplanes. Dieser umfasst ausser den drei oben genannten - 

‘ BERNOULLI noch fiinf weitere Trager dieses Namens sowie den einst beriihmten 

r Jaxos Hermann. Die Zahl der Bande witd auf 20 bis 25 geschitzt; davon 

: kommen allein auf JoHANN BERNOULLI 3-4 Bande gesammelte Schriften und 

4 6-8 Bande Kortespondenz. 


Der vorliegende erste Band erdffnet die Serie der Korrespondenzen von JOHANN 
Bernoutu. Er enthalt von Seite 95 an, in drei Gruppen A, B, C gegliedert, 
162 meist franzdsische oder lateinische Briefe aus dem Zeitraum 1690-1707. A 
* bringt alles, was noch von JOHANN Bernoutuis Briefwechsel mit seinem Bruder 
Jaxos erhalten ist. B enthilt die schon lange erwartete Korrespondenz mit dem 
Marouis DE t’Hérrrau. In der Einleitung dazu wird endgiiltig der Nachweis 
erbracht, dass die beriihmte «Analyse des infiniment petits» und andere Werke 
; pE wHérrraus fast véllig auf den Unterrichtsbriefen BERNouLLIS beruhen, was 
| ‘der Marquis bis iiber seinen Tod hinaus zu verschleiern verstanden hatte. Gruppe 
/ Cendlich ist eine Mischung nichtwissenschaftlicher Briefe von oder an rund 
i= zwanzig Personen. Im Anhang folgen noch 6 Verzeichnisse, datunter eine Biblio- 
gtaphie der Druckschriften des Marquis DE L’Hoprrat. 


i’ 2 


Zu beziehen durch Ibre Buchhandlung 


BIRKHAUSER VERLAG BASEL 


SCHUTZENMATTSTR. 31 STEINEH* to BASEL 


